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AN EXTENSION OF THE THEORY OF RESONANCES OF 
BIOLOGICAL CELLS: Til. RELATIONSHIP OF 
BREAKDOWN CURVES AND MECHANICAL Q 


EUGENE ACKERMAN 
PuHysics DEPARTMENT, PENNSYLVANIA STATE UNIVERSITY 


It has been shown that the experimental breakdown curves are much 
too sharp to be accounted for simply by a mechanical resonance of the 
biological cells. A detailed discussion of the coupling of these modes 
and their action indicate a method of estimating mechanical Q from the 
experimental curves. This value agrees with the theoretical predictions. 


In a previous discussion (Ackerman, 1954b) it was shown that 
the mechanical @ of reSonant cells should be about 2.3. On the 
other hand, the experimentally observed curves (Ackerman, 1952, 
1954a) indicate, if interpreted as amplitude vs. frequency curves, 
that the Q is as high as 6. Also a preceding paper (Ackerman, : 
1955) has shown that the coupling to these modes from a plane 
wave is completely negligible. It therefore appears worthwhile to 
examine in detail the coupling of cavitation bubbles to these modes 
and the question of the relationship of breakdown rate to 
mechanical Q. 

Coupling to Spherical Wenes: Consider first a small spherical 
source of radius, ¢, located at a distance L from the center of the 
spherical biological cell. If the center of the small source lies on 
the polar axis and if the source has a strength, S, the velocity 
potential of the source in a free field would be given by 


®, = (S/ret*r—), 
The value of ®, at the surface of the cell can be expressed as 


~ 


ee ®.(a, 9, ) = (S/L)e~*#" ) \(a/L)"P, (cos 6). 


This has no terms in P! except 7=0 and so will not excite the 


nN 
A 
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modes of interest here. If, now, the source is located in the w = y, 
plane along the line 0 = @, then we find 


®, (a, 0, w) = 
(S/L)e“*#" ) °Y “(a/L)" P, (cos 6)P!(cos 6.) cos w cos Ye. 


In particular, for the lowest resonant mode the coefficient of 
cos 2we7** is 


®,,= 98a’ L~* sin? 6, cos w, sin? 0. (1) 


If the source is a bubble driven by an applied acoustic wave, S/e 
will be of the order of magnitude of the amplitude of ®.., the ve- 
locity potential of the applied acoustic wave. In contrast to the 
plane wave case [cf. equation (6) Ackerman, 1955] if ea?L- is not 
too small, an appreciable amplitude of ®,, exists at the cell sur- 
face. In other words, a biological cell is so poorly coupled to a 
plane acoustic wave that it is only slightly effected by it. On the 
other hand, small bubbles act as acoustic transformers coupling 
efficiently to any biological cells very close to them. It is in- 
teresting to note that the magnitude of ®, at the cell surface is not 
a function of frequency. Except for resonance effects, the break- 
down rate could be interpreted as an indication of the intensity of 
cavitation. 

Active Volume. In order for breakdown to occur, the cell surface 
must. experience a certain minimum acoustic field. To drive the 


cell to the breaking point, due to vibrations in the lowest possible 
mode, 


m=2, J=2, 


the amplitude of the incident wave at the cell surface, ®,,, must 
exceed a certain critical value. Referring to equation (1) we note 
that this amplitude depends, for fixed source Strength, and for 
fixed bubble and cell radii, on the distance, L, between bubble 
and cell centers and on the angles @, and Y,. Thus, if the bubble 
center lies within a volume for which (sin? 0,cos w,/L*) exceeds 
a minimum value, B, cellular disruption becomes possible. On the 
other hand, evidence presented by Harvey et al (1946) indicates 
that the bubbles always form outside of the cell. For simplicity, 
it is assumed below that the bubble always remains spherical; i.e., 
when the bubble is against the side of the cell the center to center 
distance L is a+e, Accordingly, the bubble, effective in produc- 
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ing cellular breakdown, has its center within an active volume 
V_, defined by 


at+e<L, 
and 


—3 gin2 
L~* sin” 0, cos w, > B. 


The probability of breakdown of a particular cell will then be pro- 
portional to the number of cavities per ml, the strength of the ef- 
fective sources, and the size of V,. If the frequency of the sound 
field is varied, while holding the number of cavities and the source 
strengths constant, then the breakdown rate, K should vary as V,. 

To find the variation of V, with frequency, it is convenient first 
to find how V, varies with the parameter B. Then the variation of 
B with frequency will be discussed. Finally, these two will be 
combined, to indicate the variation predicted for V, with frequency. 

If the cell were spherical and had no preferred directions, then 
we could always choose the axes so that 


- 2 
sin“ 0, cos ¥, = 1. 


Under these conditions V, is defined by 


ate << B 


oo | 
. 


Therefore we find 
v= (4/3) 7 BO [1 — Pla +2)" 12 


However, all cells do have some preferred direction so that this 
pattern is oversimplified. 
If the polar axis, only, is predetermined, we may choose 


cos yy, = 1. 


For this case, V, is given by 
3 
V, = (8/3)7Bo1(1-Bia+ e$?)?. 


In the most general case which certainly applies to paramecia, 
all three axes will be predetermined. The active volume, V,, can 
no longer be expressed in closed form; in integral form, it is 


7 is As 3 
x ain (BL?) Des wee 


B ; 
V,=4 [ L?dL sin 640 [ be sh 
ate 1/2 
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High intensities correspond to large, V, or small B. If the in- 
tensity is sufficiently high so that 


Biase <a, 


then the active volume will vary as B™! in all three cases. If, on 
the other hand, 


Bares, 


no active volume can exist. For values slightly less than 1, all 
three cases indicate V, varying very rapidly with small changes in 
B. In the range 


0.5 <B(a+e6)* < 0.75, 


we can approximate the curve of V, plotted against B by a single 
term of the form: 
VC B 

For the perfect sphere, n= 2; for the polar axis predetermined, 
n= 3; and for both polar and & axes predetermined, n= 4. More- 
over, these formulas indicate that the shape of the curve of K 
plotted against frequency will vary with intensity. This is in ac- 
cord with the experimental observations (Ackerman, 1952), where 
it was pointed out that flatter curves were obtained with higher 
intensities and the published curves were obtained with compara- 
tively low intensities. Nonetheless, a maximum in a curve of K 
plotted against frequency must correspond to a minimum in B (or a 
maximum in B™'). The value of B in turn must be related to the 
magnitude of ®,,, the velocity potential at the actual cell surface. 
In other words, at a resonant frequency, the driving force neces- 
sary will be less. If cellular breakdown were due simply to ex- 
ceeding the tensile strength, then B should be inversely propor- 
tional to |®,, |. It seems more reasonable to suppose that, for 
these resonances, breakdown occurs due to local heating and, 
therefore, that B is inversely proportional to |®,,|? : 

To compute ®,,, it is convenient to refer to a previous paper 
(Ackerman, 1955). Equations (1), (2), and (3) in that paper can be 
used without any changes. Equation (1) of the current paper gives 


®,,- It proves desirable to lump the various constants in the 
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velocity potentials in the following way: 
®., = Cr* sin? 6, 
® , =(D/r*) sin? 6, 
®.. = Er? sin? 0, 
Substituting into the equations referred to above, we find that 


2 
lo)? = 14 


nae 1 
a? 0, a ra) ek : 


As the frequency is varied, at acoustic intensities where we can 
fit V, with a smooth curve, 


a as Bae 
we find that 


z a] - (2) 


The appearance of the power n in this expression emphasizes 
that the curve of the breakdown rate constant will be much sharper 
than a simple mechanical resonance curve, with the same quality 
factor, Q. 

Mechanical Q. As noted earlier, V, should have the same fre- 
quency dependence as the breakdown rate constant, kK. It should 
be possible using equation (2) to estimate the mechanical quality 
factor, @, of the resonance, from the experimental curves. If we 
measure the bandwidth Aw, between the two points at which K is 
decreased A times its maximum value, we find that 


1 
Q =(FA- 1)? (@,/Ao). 


The published data for P. caudatum (Ackerman, 1952), and for 
amphiuma erythrocytes (Ackerman, 1954a) have been used to com- 
pute values of Q for the cells using n= 4. These are respectively 
1.8 and 1.5. 

The actual dependence of @ on n is not very critical. If we as- 
sume n=2 the computed values of @ are raised to 2.5 and 2.1. 
Thus the breakdown curves do not make an exact determination of 
the Q of the cells possible but indicate that it is about 2 in both 
cases. This compares well with the theoretical 2.3 computed in 
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part I (Ackerman, 1954b) for viscous cells. If the choice of n= % 
were made, we would compute a value of 6 for @. The theory de- 
veloped in (Ackerman, 1954b) shows that Q is much lower than 
this. 

Breakdown vs. Intensity Curves. The theory presented in this 
paper is valid only if the applied acoustic field is approximately 
constant as the frequency is varied. If the field is increased at 
fixed frequency, the problem is more complicated since several 
effects may occur. First, the number of active centers will in- 
crease. Second, the strength of each center will be increased; that 
is, B will be decreased. Both of these will tend to increase the 
breakdown rate. On the other hand, streaming and turbulence will 
be increased. Consequently, the probability of a cell being near 
a cavity for a sufficient number of consecutive cycles will be de- 
creased. This may tend to increase B. 

Thus the treatment presented here is satisfactory to discuss the 
curves of K plotted against frequency for a constant acoustic 
field. It is also useful to discuss variations in the relative break- 
down rate of different types of cells. 

Added Comments. Other possible models would also increase 
the sharpness of the breakdown curve for a given mechanical Q. 
Thus, if the kinetic energy determines the breakdown rate we would 
have 


Ba|C|? 3. 


Likewise, it is conceivable that K might be a more complex func- 
tion of V, and that this could give a sharper curve. 

The data cannot be used to distinguish between the different 
models. However, the curves are too sharp to be accounted for by 


treating breakdown rate as linearly proportional to amplitude of 
vibration. 


SYMBOLS. USED IN II AND III 


r= radial distance 

6 = colatitude 

Ww = longitude 

® = velocity potential 

a = cell radius 

p = density 

n = viscosity 

T = interfacial tension 
subscript 7 = inside cell 


subscript 0 = outside cell (except with @ resonance) 
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subscript e = external or applied 
subscript s = scattered 
subscript a = absorption in I]; active in II 
subscript n = normal 
subscript 2 = component or constant associated with cos 2yye'®! 
og lth associated Legendre polynomial of order n 
spherical Bessel function 
spherical Neuman function 
27 times the frequency 
resonant value of @ 
velocity of sound 
o/L 
cross section 
power 
intensity 
constant associated with 9, 
constant associated with 9; 
constant associated with ®, 
sources strength of spherical sources (bubble) 
distance from center of cell to center of bubble 
radius of bubble 
cube of the inverse of the maximum effective value of L 
= volume 
number of cycles for a free oscillation to be reduced by a 
a? a 
(e] 
7, + 4 No 


sei peg | 


S) ~ a 
CoG motte SC ES ag 5 Ea F 
Houd 


tow 


factor of e~”. In part I it is shown that Q = 0.85 


for the lowest mode / = 2, m = 2. 
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MATHEMATICAL MODELS OF ENDOCRINE SYSTEMS 


LEWIS DANZIGER AND GEORGE L: ELMERGREEN 


MILWAUKEE SANITARIUM FOUNDATION, WAUWATOSA, WISCONSIN 
AND UNIVERSITY OF WISCONSIN, MILWAUKEE, WISCONSIN 


There is proposed a generalized mathematical model of endocrine 
systems, consisting of a set of differential equations which describe a 
chain of chemical reactions. The product of each reaction stimulates or 
inhibits some other reaction in the chain except possibly the last, which 
may or may not influence the system. At least one reaction must be in- 
dependent and able to proceed without stimulation or inhibition by the 
products of other reactions. 

If only two reactions of the type assumed constitute a closed chain, 
sustained periodic variations in the concentrations of the reaction prod- 
ucts cannot occur. If the chain consists of three or more reactions form- 
ing a closed loop, sustained oscillations, such as are observed in the 
menstrual cycle or in the mental disorder called periodic catatonia, can 
‘ occur under suitable conditions. In this case, the concentrations of the 
system components exhibit relaxation oscillations characterized by 
periodic degeneration of the system when an independent reaction be- 
comes completely inhibited by other reaction products. A set of con- 
ditions sufficient to produce periodicities in component concentrations 
is presented. 

Application of the model to the normally periodic system of the men- 
strual cycle and to the abnormal endocrine system which causes periodic 
catatonia is discussed. 


The endocrine system is often considered to be a mechanism 
which maintains the concentrations of certain materials within 
narrow limits. Under certain conditions, the limits are so narrow 
that the concentrations are practically constant; under other condi- 
tions, which may be normal, periodic variations of the concentra- 
tions are known to occur. 

The maintenance of the concentrations within limits is achieved 
by action of system components (hormones and activated enzymes) 
each of which stimulates or inhibits the production of other com- 
ponents. If these reactions constitute a feed-back mechanism, they 
are arranged in a closed chain, in which each component stimu- 
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lates or inhibits the production of the next component, and pos- 
sibly others in the chain. While at least one component must be 
produced independently of stimulation by another component, in- 
hibition of production of this independent component by another is 
required if the loop is to be closed. 

In such systems, a qualitative description is inadequate to pre- 
dict performance, and quantitative representation of internal rela- 
tions must be used in an analysis of the kinetics of the system. 
Quantitative descriptions of endocrine systems in the form of 
mathematical models have been proposed previously. The men- 
strual mechanism, a normally periodic system, has been investi- 
gated mathematically by Rapoport (1952); the thyroid-pituitary 
system, which may exhibit periodicities when not functioning 
properly, has been treated mathematically by Danziger and Elmer- 
green (1954, 1956). In the present paper, there is proposed a 
generalized mathematical model and there is formulated a set of 
conditions sufficient to yield periodicities in the component 
concentrations. 

Consider a system of n components whose concentrations in the 
body are functions of time and are denoted by 2,1, 2,..., #,. The 
Kinetics of the system may be described approximately by the set 
of first-order differential equations 


de /dt+h ;,=Q;,  (¢=1, 2,3,...,2) (1) 


wherein 4, is the loss rate per unit concentration of the component 
z,, and Q; is the production (activation) rate of the hormone (enzyme) 
a; The production (activation) rate is, in the first approximation, 
a linear function of the system components and may be represented 
by 


Q; = Aig + a A;;2;, (2) 


j Poa 


wherein A; is zero or a positive constant denoting independent 
production (activation) of 2,;; the A ij are sensitivity constants and 
may be zero for no effect, positive for production stimulation (acti- 
vation), or negative for production inhibition (deactivation). The 
constant A;; is considered to be zero in all cases as auto-stimulation 
or inhibition is not likely. 

The function Q; as given by (2) may take on negative values for 
some positive sets of values of the a; If 2, is a hormone, a nega- 
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tive production rate is not possible and the restriction Q@.> 0 must 
be imposed; from (1) this also limits 2; to non-negative values. Me 
a; 1S an enzyme, a negative Q; represents deactivation and is 
permissible; in this case the restriction is simply z,>0 as nega- 
tive concentrations of activated enzymes are not physically possible. 

With these restrictions, the mathematical model may be formulated 
completely as 


dx ;/dt + d;x; = Ajo + a Ajj; = Q; 
j=1 


w,20, @Q,;>0 if; is a hormone (3) 


ee Soh ees 


Rapoport (oc. cit.) includes, in his model, the possibility of 
negative production rates of system hormones by assuming that one 
hormone can enhance the dissipation or breakdown of another hor- 
mone in the blood stream. This assumption leads to linear equa- 
tions but, as Rapoport points out, is not in accord with observa- 
tion. The model (3) represents inhibition action wherein com- 
ponents act inhibitively on the cells which produce hormones; this 
assumption precludes hormone production at negative rates. Ac- 
cordingly, the model presented here will be characterized by the 
degeneration of one or more hormone equations when the associ- 
ated Q; are driven to zero. Degeneration of one equation, in ef- 
fect, isolates the mechanism which produces the hormone from the 
influence of the rest of the system. In the interval during which 
(2) yields a negative value, the degenerate equation is a first- 
order, homogeneous, linear differential equation whose solution is 
an exponential decay of z; with a decrement A,. 

The degeneration of the mathematical model produces a system 
which is not linear in the unrestricted sense. Instead it is piece- 
wise linear comprising a set of linear equations which may pos- 
sibly be unstable in the normal state and a different and stable set 
of linear equations in the degenerate state. 

The concept of a ‘‘push-pull’’ endocrine system has been applied 
to the mechanism of the pituitary and thyroid glands in the regula- 
tion of the metabolic rate and to the pituitary and ovary in the con- 
trol of the menstrual cycle. The term ‘‘push-pull’’ indicates stimu- 
lation or push by the first gland in the closed loop of production by 
the second, and inhibition or pull by the second gland of produc- 
tion by the first. Clearly, the ‘‘push-pull’’ system exhibits the 


12 LEWIS DANZIGER AND GEORGE L. ELMERGREEN 


phenomenon of feedback. The simplest ‘*push-pull’? system can 
be represented in the form (3) by the following second-order set in 
which the a;; are positive constants: 


dx, /dt +d 121 = G19 — ort, = Qi, Q,>0 


(4) 


The set (4) reflects the following assumptions: 

(1) only 2, is produced independently. 

(2) w, is a hormone and cannot be produced at a negative rate. 
The system is degenerate for x > @49/@12- 

(3) x, stimulates the production of the hormone 2p. 

(4) az. inhibits the production of 2,. 

(5) the components dissipate in proportion to their respective 
concentrations. 

As the system (4) exhibits a degeneration point at x = @j9/@4o, 
different solutions exist for the normal system with ag < a, )/d@12 
and for the degenerate system with zg > @j)/a@,,. These two solu- 
tions can be obtained by standard methods and pieced together at 
the degeneration point, x, = @49/d@,,. Each solution has the form 


2 
a, = 6; + 3 C;; exp (7,2), (4 = 1, 2) (5) 
j= 
where the 0; are the steady states, the c,; are integration con- 
stants, and the r, are the roots of the characteristic equation 
[multiple roots being excepted if the form (5) applies]. 

The possibility of sustained oscillations of the component con- 
centrations in a “‘push-pull’’ system can be investigated by ex- 
amining solutions of typical systems. The second-order system (4) 
for the normal condition, z, < a@4)/aj2, has steady states which are 


by = AyAg/(A Ay + @12%a1) 5 


(6) 
bg = Gyo4a1/(AyAg + AyQ4q}) 5 
while the 7, are the roots of the characteristic equation 
(D + 24) G19 0 ~ 
= Go% (D + dg) a 


For the degenerate condition, x, > @19/@,., the system differs only 
in that @, = 0. Equivalently, set a1) = a1) = 0 in (6) and (7) yield- 
ing, respectively, 6, = b2=0 and r, = -A,, 7 = ~A,. In the 
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normal state, the roots of (7) have negative real parts for all posi- 
tive values of the coefficients, and solutions for x < a,)/a,) must 
either be aperiodic or be a damped periodic variation about the 
steady states of (6). The steady-state concentration, b,, from (6) 
is less than the value of xz, at which degeneration occurs; the 
steady states for the normal condition are, therefore, stable equi- 
librium points. The degenerate state, x > a@,)/@,2, would obtain 
only during an initial transient and then only if an initial condition 
to(t = 0) > aio/ayq exists. The degenerate state, characterized by 
zero steady states and negative, real r,;, would exist only until z, 
fell below the degeneration point after which the normal state 
would exist for all time. Clearly, the second-order system (4) 
could not exhibit sustained oscillations of the component 
concentrations. 

If, in the system (4), the normal state had been unstable so that 
the component concentrations would tend to increase without limit, 
the system would degenerate regardless of initial conditions. Re- 
laxation oscillations would then exist as the system would alter- 
nate between the unstable normal state and the stable degenerate 
state. Although this is not possible in the system (4), a physically 
realizable set of conditions which would produce sustained oscilla- 
tions in a similar manner is proposed for ‘‘push-pull’’ systems of 
order greater than two. These conditions are: 


I. at least one of the 7; for the normal system has a positive 
real part; 

II. all the r; for the degenerate system have negative real parts; 

III. the steady-state concentrations for the degenerate system 
are less than the concentrations at which degeneration occurs. 
If these conditions obtain in a ‘‘push-pull’’ system, relaxation 
oscillations will occur with the following sequence: 


(A) with the system in the normal state, the roots, r;, which 
have positive real parts would cause the system to be unstable, 
and all concentrations would ultimately rise. This unstable state 


would persist until the system degenerates. 
(B) when degeneration occurs, the system enters a stable mode 


of variation for which all the 7; have negative real parts and would 
tend toward steady states which are lower than the concentrations 
at which degeneration occurs. The concentrations in this stable 
mode would ultimately fall until the unstable normal state would 


again apply. 
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The processes A and B would repeat cyclically, and the system 
would be governed alternately by the unstable and stable solutions 
of the normal and degenerate state equations. The resulting os- 
cillations would be established for any set of initial conditions. 

As_was shown, the second-order system (4) cannot meet condi- 
tion I; further no second-order system of the form (3) can meet the 
requirement of an unstable normal state as the 7; must have nega- 
tive real parts if the loss constants, A;, are positive. 

We will examine next the simplest third-order ‘‘push-pull’’ system 
which has the form 


adr ,/dt +12, = Ayo — Ayg%3 = Q,, @,>0 
dx, /dt + Ag te = 221%1, (8) 
If the a;; are positive constants, the set (8) reflects the following 
assumptions: 
(1) only x, is produced independently. 
(2) w, is a hormone and cannot be produced at a negative rate. 
The system is degenerate for 2, > @49/a13- 
(3) x, stimulates the production of 22. 
(4) x2 stimulates the production of z,. 
(5) xz inhibits the production of z,. 


(6) the components dissipate in proportion to their respective 
concentrations. 


For the set (8), the normal system, with x3 < a,)/a,,, has the 
characteristic equation: 


(D + X,4) 0 Qis 
—Qo1 (D + Ag) 0 =0. (9) 
0 — 239 (D + Xs) 
The first of the proposed conditions for sustained oscillations re- 


quires that (9) have at least one root with a positive real part. If 
(9) is expanded as 


A,D® + A,D? + A,D + A, =0, (10) 


the stability criterion of Routh, see Evans (1954), applies and 


Shows that condition I is satisfied if the physically realizable 
inequality 


AyAs > 4,45 (11) 
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obtains. The second and third conditions apply to the degenerate 
state which differs from the normal state only in that @,=0. 
Equivalently, set @,) = a@,, =0 in (8) and (9) yielding r, = -r,, 
Tg = ~Ag, 73 = —Ag. As the r,; for the degenerate state are all 
negative and real, condition II is satisfied for all positive values 
of the loss constants. Condition III is also satisfied for all posi- 
tive values of the coefficients as with Q@, =0, 2, decays expo- 
nentially toward a zero steady state, and z, and z, follow. If the 
inequality (11) is satisfied, the system (8) will meet all three of 
the proposed conditions and relaxation oscillations will occur. 

Illustratively, consider the hypothetical but possible third-order 
system: 


dx,/dt+0.52. =2,, (12) 
dx,/dt+1.5a, = Ta. 


If the determinant (9) is expanded with the coefficients of (12), 
the characteristic equation 


D3 + 3D? + 2.75D + 63.75 =0 (13) 


is obtained. Applying the inequality (11) to the coefficients of 
(13) yields 63.75 > 8.25, and condition I is satisfied by the system 
(12). The existence of a root with a positive real part may be 
verified by factoring (13) as 


(D + 5)[(D — 1)? + 11.75] =0. (14) 


The 7, from (14) are —5 and (1+ 73.43), and the 2, solution for the 
normal state may be written in the form of (5) as 


= 13. —-13.43t 
La = 0.1098 + C31 € wie | e’ feet aa + Cag € ; ie (15) 


where zx, < 1/9. 
For the degenerate state we have 


- —t 70.6.5 2% lb t 9, 16 
By = Cai | + Chae + Cya€ : el (16) 


The solutions (15) and (16) would apply alternately. With 2, < 1/9, 
the pair of complex roots of (14) yield a periodic term with expo- 
nentially increasing amplitude and a3, from (15), would ultimately 
reach the concentration of 1/9 unit. At this point, dx,/dt is posi- 
tive and z, would continue to rise to a maximum value which is 
greater than 1/9. During the time that a, > 1797 €16) applies and 
a, would ultimately fall toward the zero steady state until the 
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transition point, 2, = 1/9, is reached. At this point, dz;/dt is 
negative and x, would continue to fall to a minimum value which 
is less than 1/9. During the time that x, < 1/9, (15) again applies 
and the cycle repeats. A particular solution for the sustained 
oscillation would require that the integration constants in (15) and 
(16) be evaluated at the transition point, z, = 1/9. 

The third-order system (8) has been employed by Danziger and 
Elmergreen (loc. cit.) to describe the thyroid-pituitary system 
wherein x,is the concentration of pituitary hormone, thyrotropin; zx 
is the concentration of an activated enzyme within the thyroid 
gland; and 23 is the concentration of thyroid hormone. The peri- 
odic variation of metabolic rate observed in periodic catatonia, a 
mental disorder, was explained by the relaxation oscillation theory 
presented here. Verification of the form of the solution was ob- 
tained using analog computer methods wherein adjustment of 
system parameters to obtain periodicities was possible; the re- 
quirement imposed by condition I and the inequality (II) was 
confirmed. ; 

Rapoport (loc. cit.) gives two models for the pituitary-ovary sys- 
tem which controls the menstrual cycle. If negative hormone pro- 
duction rates are not permitted, the following fourth-order set will 
include both of these models: 


dx, /dt + A121 = Ayo — Aygg — Aye, = Qi, Q,>0 
dx, /dt + rote = Ag Ly ) 
dz, /dt + Ax, = Ago 


dx,/dt + Ag@, = G43%3 . 


(17) 


Bo 


If the a;; are positive constants, the set (17) reflects the following 
assumptions: 
(1) only a, (the pituitary hormone, prolan A) is produced 
independently. 
(2) #, cannot be produced at a negative rate. The system is 
degenerate for @y9% + @142%4 > Qo, 
(3) a, stimulates the production of x, (the estrogenic hormone 
of the ovary). 
(4) vw. stimulates the production of a, (the pituitary hormone, | 
prolan B), 
(5) a, stimulates the production of a, (the luteal hormone of 
the ovary). 
(6) both x and x, inhibit the production of iis 
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(7) the components dissipate in proportion to their respective 
concentrations. 
For the set (17), the normal system, with Ayo%y + Ay4%4 < Ay, 
has the characteristic equation 


(D + X,) Gis 0 Q14 
ao (D + AQ) 0 0 So" 7is) 
0 — Ago (D + As) 0 
0 0 — A483 (D + X4) 
If (18) is expanded as 
A)D* + A,D? + A,D? + A,D + A, =0, (19) 


Routh’s criterion (loc. cit.) yields the physically realizable 
inequality 


Jia A Asi Ah -A.. (20) 


Condition I is satisfied in this fourth-order system if (20) obtains; 
conditions IJ and III are satisfied for any set of positive coefficient 
values by the same arguments that applied to the set (8). In addi- 
tion to being of higher order, the set (17) differs from (8) in that 
degeneration occurs along the line a,5%2 + @14%4 = @yq instead of 
at a point. This occurs because (17) represents a feed-back sys- 
tem with two loops instead of one as in (8). Application of Routh’s 
criterion to (17) yields two conditions for instability for the normal 
system: however, only the condition of (20) is physically realizable 
with positive a;;. 

We believe the system (17) is a good mathematical approxima- 
tion for the control system of the menstrual cycle. Verification 
would require the determination, by experiment, of the system co- 
efficients. A solution of (17) with numerical coefficients would 
then yield the time variation of component concentrations and the 
period of oscillation. 


General Remarks 


The mathematical model (3) is an improvement on the model of 
Rapoport (loc. cit.) in that it describes the steady-state behavior 
of endocrine systems exhibiting periodicities in component con- 
centrations. The essential difference between this and earlier 
models is the exclusion of negative hormone production rates and 
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negative component concentrations. Degenerate equations result 
producing a limiting action which ensures conditional stability so 
that component concentrations are bounded even if the system is 
unstable in the normal state. Although the model becomes non- 
linear if the system degenerates, the equations yield to simple 
analysis from which an existence criterion for periodicities is 
established. If time solutions are required for periodic systems 
for which numerical values of coefficients are known, analytical 
methods will yield only approximations. The use of differential 
analyzers, analog computers, or other machine methods of solution 
would be desirable. 

Work in progress suggests that the model (3) will also describe, 
to a first approximation, the mechanism which regulates the blood 
sugar concentration and the mechanism of the pituitary-adrenal 
system. 


The authors are much indebted to Drs. Anatol Rapoport and 
Ernesto Trucco for valuable criticism of the manuscript. 

This work was aided by a grant from the University of Wisconsin 
Alumni Research Foundation. 


LITERATURE 


Danziger, L. and G. L. Elmergreen, 1954. ‘‘Mathematical Theory of 

Periodic Catatonia.’? Bull. Math. Biophysics, 16, 15-21. 
and - 1956. ‘*The Thyroid-Pituitary Homeostatic Mech- 

anism.’’ Bull. Math. Biophysics, 18, 1-13. 

Evans, W. R. 1954. Control System Dynamics. New York: McGraw-Hill 
Book Co., Inc. 

Rapoport, A. 1952, ‘*Periodicities of Open Linear Systems with Posi- 
tive Steady States.’’ Bull. Math. Biophysics, 14, 171-83. 


RECEIVED 8=2=-56 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 19, 1957 


TOPOLOGICAL BIOLOGY: A NOTE ON RASHEVSKY’S 
TRANSFORMATION T. 


ERNESTO TRUCCO 


MENTAL HEALTH RESEARCH INSTITUTE 
UNIVERSITY OF MICHIGAN 


In the bio-topological transformation between graphs denoted by «Toe Xx) 
N. Rashevsky (Bull. Math. Biophysics, 18, 173-88, 1956) considers the 
number of fundamental sets which (a) have only one specialized point as 
source (and no other sources), (b) have no points in common (are ‘‘dis- 
joined’’); he proves that this number is an invariant of the transformation. 
In this note we show that Rashevsky’s Theorem can be extended as fol- 
lows: The number of fundamental sets of the first category is an invariant 
of the transformation. We must, however, count the subsidiary points of 
the transformed graph as specialized points. We recall that fundamental 
sets of the first category are those whose sources consist of specialized 
points only (Trucco, Bull. Math. Biophysics, 18, 65-85, 1956). But in 
this modified version of the Theorem the fundamental sets may have more 
than one source and need not be disjoined. 


In studying some properties of his bio-topological transformation 
T, N. Rashevsky (1956, pp. 185-87) has proved the following theo- 
rem: The number of fundamental sets (FS’s) which (a) have only one 
specialized point (and no others) as source and (b) are disjoined is 
an invariant of the transformation; i.e., it is the same for the trans- 
formed graph 7(P) as for the primordial graph P. In this note we 
show that Rashevsky’s theorem can be generalized in two ways, 
namely: (1) the FS’s in question may have more than one source, 
provided only all these sources are specialized points (or subsidi- 
ary points in the transformed graph), (2) the FS’s need not be dis- 
joined, 

In the transformed graph, 7, we distinguish between specialized, 
subsidiary, and nonspecialized or residual points. If we use the 
term ‘‘nonresidual’’ to denote both specialized and subsidiary 
points, we may classify the FS’s of T as follows: 

a) FS’s of the first category; their sources consist only of non- 
residual points, — 
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b) FS’s of the second category; all their sources are residual 
points, 

c) FS’s of the third category; their sources include both non- 
residual and nonspecialized points, 

The same classification applies to the primordial graph P, in 
which, however, ‘‘nonresidual’’ and ‘‘specialized’’ mean the same 
thing, since there are no subsidiary points. In a previous paper 
(Trucco, 1956, referred to as I in the following) we called p’4, p’4 
and p’’’ the number of FS’s of the first, second, and third category 
in P, 

With the above terminology our theorem may now be stated as 
follows: The number of FS’s of the first category is an invariant of 
the transformation T. 

As in I we denote by P, and P, the graphs obtained after the 
steps T, and 7, of the transformation, Also, we call s‘P) the 
specialized points of the primordial graph P, and s, (or s,.) the 
corresponding points of P,; the subsidiary points attached to s_. 
are called s,,(v 21). The index 7 runs from 1 to n, and the in- 


dex v from 0 to /,; the J,’s are non-negative integers such that 
n 


>» 1, = n(m — 1). 
j= 

The correspondence between the point-bases of P and P, was 
made clear in I, Furthermore, each non-subsidiary point of P, is 
a source if and only if it corresponds to a source of P,. However, 
some of the subsidiary points in P; may also become sources. In 
particular, if the specialized point s; is a source of the third cate- 
gory in P, all its subsidiary points s;, are also sources, whereas 
if s; is not a source none of its s,,’s can be a source (this follows 
easily from the rules of the transformation 7 and from Theorems 1 
and 2 of I). 

There is thus a one-to-one correspondence between the FS’s of 
the first category in P and in P,. The same is true for the FS’s of 
the third category; the number of FS’s of the second category, on 
the other hand, is multiplied by m in going from P to P,, so that the 
point-base number of P, has the value p’-+ mp’ + p*’’. 

We still have to consider the last step of the transformation T, 
1.0.5 the addition of n(m — 1) residual graphs Rye one such graph 
being attached to each of the subsidiary points. 

First of all, since the residual graphs consist of nonspecialized 
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points only, their presence cannot add any specialized sources to 
the final graph, and therefore the number of FS’s of the first cate- 
gory, p’, will not increase. 

Next, let I] be a FS of the first category in P, and B the point- 
set formed by its sources. All points of B are nonresidual, and 
any two of them are connected by at least one way. Denote by 
B, » the point-bases of the residual graphs Re The FS II will 
Boase to exist and become part of some other FS of the second or 
third Sites ea if and only if there is a way going from at least one 
of the B, y Sy Say is Pe to B (i.e., from a point of B;, , to a point of B). 

Since ‘the connection between B%; and P, is made through a sub- 
sidiary poroly say s%,, this implies that: (a) there is a way from 
By to s%,, (b) the subsidiary point s‘, is a source of II (I, Theo- 
rem 2), 

This, in turn, means that the specialized point s;-must belong to 
II and also be a source of II. But then condition (a) is not satis- 
fied (I, end of page 80), and hence the number of FS’s of the first 
category remains constant throughout the transformation 7, which 
is what the theorem states. 

Moreover, let II,‘”?, II,‘”? be two FS’s of the first category in 
P, Ul,, U1, their corresponding FS’s in 7. Then it is easy to see 
that Il, and I, will be disjoined if and only if the same holds for 
I, (P) ae ibe CP), 

ee mors: if a FS II‘? of the first category has only one spe- 
cialized point (say s,‘”) as source, the same is true for the cor- 
responding FS II in the transformed graph: none of the subsidiary 
points s,, attached to s; can become a source (this does not hold 
for some ‘other types of transformations considered by N. Rashevsky). 

In this way we obtain Rashevsky’s special case of the theorem. 
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In this paper the characteristics of a membrane having properties simi- 
lar to those of frog skin are considered in relation to renal function, Such 
a membrane, as discussed by H. H. Ussing and K. Zerahn (Acta, Physiol. 
Seand., 23, 110, 1951), is capable of acting as a sodium pump actively 
transferring sodium ions through the membrane. Water transferis regarded 
as dependent upon osmotic and hydrostatic pressures and on the permea- 
bility of the membrane to water. 

It is shown that for a membrane having these properties the rate of 
transfer of sodium is a function of the rate of transfer of water and vice 
versa, In this case a change in, say, the flux of water across the mem- 
brane may lead to a similarly directed change in the sodium flux across 
the membrane. If the E.M.F. developed by the sodium pump is controlled 
by a sodium retaining hormone (SRH.) and the permeability of the mem- 
brane to water by antidiuretic hormone (ADH.), it follows that there is an 
interaction between the two hormones which is a property of the membrane 
upon which they act. 

It is suggested that the mechanism of action of the renal tubules with 
regard to sodium and water excretion may be not entirely dissimilar to 
that discussed for frog skin, and on this supposition the characteristics 
of the renal response to certain physiological stimuli are deduced and 
found to be in excellent agreement with observation. In particular, the 
excretion of a urine initially hypotonic in response to an isotonic saline 
infusion is explained as is also the inability of, physiological doses of, 
ADH to inhibit isotonic or hypertonic saline diuresis. 


I. Introduction. 


During investigations of the renal excretion of water and sodium 
by rats which had been infused with saline it became necessary to 
consider possible modifications of renal tubular water reabsorption 


*Dr. Cole is in the Department of Medicine, and Mr. Meredith is in the 
Department of Electrical Engineering. 
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by that of sodium, and vice versa, Although in the following paper 
we will consider a model system, rather than renal function as 
such, it is not unreasonable that the behavior of such a model may 
usefully be considered when trying to account for the vagaries of 
sodium and water excretion. 

H. H. Ussing and K. Zerahn (1951) conceived of sodium trans- 
port across frog skin as being related to the activity of a ‘*sodium 
pump’’ which was regarded as an equivalent E.M.F. transferring Na 
ions through the membrane in the face of the membrane resistance 
and the electrical potential difference across it, when the two sides 
of the membrane formed parts of a concentration cell. This sodium 
pump was an active, energy consuming process located in the mem- 
brane itself. Water transfer was regarded as dependent on osmotic 
and hydrostatic pressures across the membrane and on the permea- 
bility of the membrane to water. In the case of the kidney there is 
no evidence of tubular secretion of sodium or water, and we may 
consider the mechanism concerned with the net transfer of these 
substances to act unidirectionally with regard to their reabsorption. 
For the present, we have disregarded the direct effects of water 
flow on Na transfer (Teorell, 1953) and propose a model in which 
the E.M.F. developed by the sodium pump is controlled by a sodium 
retaining hormone (SRH) and where the membrane permeability to 
water is controlled by antidiuretic hormone (ADH) (Sawyer, 1951; 
Keofoed-Johnsen and Ussing, 1953). We propose to show that the 
nature of the membrane can account for some apparent interaction 
between SRH and ADH, so that this interaction will depend rather 
on the target-organ of the hormones than on any direct interaction 
between them. We do not underestimate the importance of direct 
action of one hormone on another; e.g., an effect of ADH on SRH 
production, but we do intend to consider first some possible con- 
sequences of interactions dependent only on the properties of the 
target organ. 


Il. Behavior of a simple artificial system. 


We shall commence by considering some features of a simple 
model and discussing its equilibrium states. Certain properties of 
this system will be described and may be used in the consideration 
of more complex, and more nearly physiological, systems. 

Consider a membrane separating two compartments, A and B, 
each containing aqueous solutions of sodium chloride. 
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Glomerular 
filtrote 


tubular tia 6 steqdy state peer nae 


4 Cus Compt. B -E.C.F 
ae Compt. A- 


Urine 
FIGURE la. Representing Na and water fluxes between two compart- 
ments, A and B, separated by a membrane whose permeabilities to Na and 
water are Cnq and Cy respectively. S, is the quantity of sodium in 4; 
V4 the volume of water in A. 
FIGURE ib. Representing a simplified model of a renal tubule whose 
reabsorptive properties are uniform over the whole of its length. 


Symbols: 


V a(2) = volume of water in compartment A at time ¢ 
Va(2) =z ”? ” ” ” ” B at time ¢t 
S(t) = quantity of Na’ in A at time ¢ 
Sg(t) = quantity of Na’ in B at time ¢ 


Cnq = netmembrane permeability to Na’ (in the direction A—> B) 
Cw =net membrane permeability to water (in the direction 
A — B) 
ina = net Na’ flux, A—» B 
zt» = net water flux, A— B 


The symbol Na’ represents sodium ions. © 
It is also supposed that V4 is much smaller than Vg, so that the 
error in taking V4 + Vg = Ve may be neglected. 


A. Where there is no active transfer of either Na or water. 
At time ¢, the force tending to transfer Na* from A to B depends 
on the logarithm of the concentration ratio and is expressed by: 
: Sa(d)/V a(t) 
Sp(t)/Ve(t) 


Sa(d) Sp (t) 
i.e., K {108 5" (2) — log 87 ee (1) 


where K is a constant. 
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Similarly, with equal hydrostatic pressures on the two compart 
ments, the force tending to transfer water from A to B may be ex- 


pressed as: 
Sp(?) Sa(?) 
K {i Va(2) — log Be (2) 


The flux for the two substances is given by the product of force 
and permeability, hence: 


Sa(t) SB(¢) 
tNa ~K {og Va(t) — log ae (3) 


and 


Sp(t) Sa(d) 
ty le {oe V(t) - log aed Cas (4) 


assuming that the membrane permeability is not a function of flux. 
Now, 


; adSa4 dSpg 
eeteeds Sg 
and 6 ge 
at at 
And since Vg >> Va, (3) and (4) become at time ¢: 
Sa(t) Sp dS 
K « log —— - log —}Cn, = —-— 
{ Viléin a ae at’ (5) 
and 


S 
K {tog ane log 
B 


V, Va(t) 


Sa(t) av 
af bo. =~. (6) 


Because the transfer is unidirectional we add the condition that 


dS dV 
oe <0 and ors < 0. (7) 


From an examination of equations (5) and (6) it is clear that an 


nie re Ss 

equilibrium condition exists where ee a ‘whatever the initial 
A 

concentration in A, Equilibrium is reached when there is no con- 


centration difference across the membrane. 
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B. Where there is active Na transfer. 


We now proceed to consider the modifications of the system just 
described when the membrane itself acts as a sodium transporting 
system (‘‘sodium pump’’) in the direction A — B (cf. Ussing and 
Zerahn, loc. cit.), while net water transfer remains a permeability 
effect as before. As in the work of Ussing and Zerahn (loc. cit.) 
the action of the ‘‘sodium pump’’ is expressed in terms of a poten- 
tial, Ena. 

Equation (5) now becomes: 


dS 4 Sa(d) Sp 
eee net KA) Steet One, 8 
[En : {ler EVs [aloe 2 


while (6) remains unaltered. 


dS 
It is clear from equations (6) and (8) that with Ey, > 0, a and 


aVa 
— cannot simultaneously be zero and that, therefore, no equilib- 


rium condition exists comparable with that of the previous case. 

However, a steady-state condition may be shown to exist such 
that avs and pee take non zero values which keep Sal?) 

dt dt Va(t) 

fore sodium and water flux from A to B) constant. There is, then, 
no equilibrium but there is a steady-state condition. Thus there 
will be complete transference from A — B unless there is an op- 
posing force such as hydrostatic pressure. 

In the steady-state condition we have: 


(and there- 


Sa(t dS dVA 
A(é) = constant, er fe constant, and —— = constant. 
Va(e) dt dt 
Sa(t 
Since ) is independent of time (= «), 
A 


Sa(t) = aVa(?). 


Differentiating with respect to time: 


dS, e dV4 
ES see ae 
dS 
S dt 
Se (9) 
V4 d~ 
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aV, 
Substituting for ae from equation (6) and (9) in (8): 


S(t) SB aS a 
EK =lor=— de Le eI 
CNna | NatK {oe Va(d) og Vp 7 


Sp aA 
= a-K jlog — - log ——} Cu; 

2 Vz = Va(2) 
Sa(té S EN Cre 

log 4 og 
Va(t) Ve K(Cna + XCw) 

Sa(t) Sp. e 
From (10) it follows that ——— < — in the steady-state condition; 
Va(t) Vg 


that is, the compartment A is hypotonic with respect to the com- 
partment B, and its actual concentration will depend upon Ewa, 
Cnq and Cy. For example, if Ena is increased the steady state 
would alter so that A became more hypotonic with respect to B 


2-Ot Hypertonic 


Hypotonic 


ea) 
O25 
fe, 
o Gs, P=O 
i 
Cc 
g Ol 
2 
ic 
e 0:08 E,=7OmvV 
3 
O 
Cit terme Sf oes 
0:06 
OoIc,. O1CA, tOC\, 1OCyg 
Permeability to Water 
FIGURE 2, 


The effect of Eyq upon steady-state concentration ratio 
( 
between A and B, me plotted as a function of C, with P = 0. Note that 


B 
compartment A is hypotonic to compartment B throughout, except where 


Ena =0, and that the degree of hypotonicity increases with increasing 
values of Hyg. For details of parameters see text. 
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(Fig. 2) and, from (6) supposing Cy, and Cy unaltered, the flux of 
water from A to B would also increase, 
C. Where there is a force favoring water transfer. 


Suppose that there is a force, P, tending to transfer water from 
A to B (e.g., hydrostatic pressure, or osmotic pressure due to a 
solute other than NaCl), equation (6) is modified: 


aV4 Sp Sa(?) 
Se + Kh log —— = | 
dt k | SS a wi mal Oe 
ae d 
Again, no equilibrium condition exists since “ and <4 cannot 
t 


simultaneously be equal to zero. Proceeding as before to evaluate 


Sa(t) . 
in the steady state with Ey, =0; using equations (5), 


Va(d) 
(10), and (11) we have: 
dS S 4(2) SB 
Ber aK op = log —} Owe 
dt | Vag Vi (12) 
Sp S(t) 
=a |P+K<+log —-lo (ORS 
| Ve VAC?) 
S4(t) Sp aP Co 
lo OS ar ae, 13 
© V A(t) e Ve K(Cna + X0w) ue 
S(t) : Sp 


> — in the steady-state condition 
Va(t)” Ve 

so that compartment A is hypertonic with respect to compartment 
B, at least with regard to Na’ concentration, (Fig. 3). 


From (13) it follows that 


D. Where there are forces favoring transport of both Na and water. 

In the presence of both these effects together, that is, where 
both P and Ey, > 0, the steady-state equation is: 
S4(o) Sp &PCyw-Ewna Cna 
Vie olan (Ole abe) 
At the outset we supposed that the variable properties of the 


membrane were Ey, and Cy. It is evident from equation (14) that 
in the steady-state condition the contents of A may be either hypo- 


(14) 


log 
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A:B 


ratio 


Concentration 


Permeability to Water 
FIGURE 8. The effect of P upon the steady-state concentration ratio 


between A and B, =e plotted as a function of Cy with Eng =0. For 


details of parameters see text. 


E.= 2OmV 


A:B 
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FIGURE 4a, The effect of P upon the steady-state concentration ratio 


between A and 8, —— , plotted as a function of Cy with Evg = 20 mV. 
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FIGURE 4b, The effect of P upon the steady-state concentration ratio 
Ca 
between A and B, cea , plotted as a function of Cy with Eyg = 50 mV. 


For details of parameters in these figures, see text. 


tonic or hypertonic with respect to B and that this can depend on 
the values taken by Ewa, and Cy (cf. Figs. 4a and 4b). 
Sa(é) SB 


Substituting for {log — log — 
Va(t) 8 Vp 


} in equations (8) and (11) 


from (14), we obtain: 


adS4 aPCyw-EnaCna 
te Ae Se ee eee ee ee 15 
= {ene Ggerk Os hone, (15) 
and 
w-Ena C 
s dV = a GPG Na UNa pe (16) 
dt Cnat+ACy 
; dS 4 Ati 
These equations give the values of pe and rar in the steady 
state. 
From (15) and (16) it is evident that a change of Eq alters the 
d V, 
values of both “2 and — similarly a change of C,, alters both 
= and ee There is thus an interaction between Fy, and water 


32 D. F. COLE AND J. F. MEREDITH 


transfer, and between C,, and sodium transfer. If we consider the 
effect of varying Eva and C, on the steady-state condition a series 
of ‘‘trajectories’’ may be drawn from which this phenomenon of in- 
teraction may be seen more clearly (Fig. 5). 


SODIUM REABSORPTION C Increasing excretion) —> 


dS, 
dt C 
755) 20 15 IO 5 @) C= Neve) 
WATER 


Me REABSORPTION 


(Increasing excretion) 


© —> E,= 20 mV 
@ — £=50 mV 
A —» £,= 70 mV 


dS 4 dV, 
FIGURE 5. Graphs of ra against 


with Ey, a continuous and C, 


a discrete parameter. Each trajectory is drawn for a constant value of 
Cy. Symbols: Eng = 20 mV; Eng = 50 mV; Eng = 70 mV. If these curves 
represent conditions in a tubule (cf. Fig. 1b), then decreasing reabsorp- 
tion corresponds to increasing excretion, subject to the limitations dis- 
cussed in the text. For details of parameters see text. 


Values of constants and parameters 


a = nee unit concentration being taken as i, 
Va B 
Cw is a parameter: with a NaCl ratio of 1:10 across frog skin, 
the water flux is about 10n1/cem?/hr. (Koefoed-Johnsen and 
Ussing, loc. cit.; Capraro and Garampi, 1956). From these 
figures, C,, may be expressed as a membrane conductivity 
of 2.7 x 10-° ml/min/cm?/mV. In the presence of ADH, the 
conductivity C,, will be more than doubled. 


Ena is a parameter: Ussing and Zerahn (loc. cit.) give values of 
about 70 mV for frog skin. 
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: RT 
K is a constant and = is (where n = 1) and is taken as 26.6 
n 


mY. 
Cna is the same as kya used by Ussing and Zerahn (loc. cit.) 
and is equal to 3.7 x 10-® mEq/min/cm?/mV. 
P is a parameter of the system given, for convenience, in 
millivolts, 


Ill. Discussion of applications to renal tubular function. 


Suppose we consider a simplified model of the renal tubule in 
which the reabsorptive properties are uniform over the whole of its 
length and follow the progress of an element of solution entering 
at the glomerulus and flowing downwards. A part of this element 
of solution is reabsorbed into the extracellular fluid and vascular 
system (corresponding to compartment B), while the remainder en- 
ters the ureter as urine. The element of fluid will have been in 
contact with the tubule wall, whose properties are regarded as 
analogous to the membrane in the model, for a time 7, equal to the 
time taken for the element of fluid to flow down the length of the 
tubule. This is similar to the situation in the model where the 
fluid remaining in the compartment A (i.e., the tubule lumen) at the 
time ¢ = T is termed urine. E. Bojeson’s experiments (1954) sug- 
gest that in the kidney the value of 7 is not much less than 60 
sec. At present it is not proposed to discuss the morphology of 
the renal tubule or its reabsorptive processes but it is possible to 
consider the case where a ‘‘sodium-retaining hormone’’ (SRH) con- 
trols active Na transfer by its action on Ewyq and where tubular 
permeability to water (C,) is increased by the action of ADH. 
This latter view is not dissimilar to that proposed by W. H. Sawyer 
(1956). 

The present approach will imply that the hormones concerned 
not only cause reabsorption of specific ions or molecules but also 
alter the steady-state reabsorption condition of both Na and water, 
and we propose to consider how such a model would behave in a 
number of physiological conditions. 

_ In this discussion it will be assumed that: 
(i) the time taken for an element of fluid to traverse the length 
of a tubule is constant and equal to 7, 
(ii) the steady state of reabsorption is achieved in a time much 
less than T. 
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Na REABSORPTION —<— 


WATER REABSORPTION 


FIGURE 6. (Reduction of sodium reabsorption) For any trajectory a 
reduced reabsorption of Na’, due to a fall of E Na, leads to increased water 
excretion. Arrows on the trajectories represent decreases of Ena. 


: : ; dS, dV4 
In this case the axes of Figure 6, i.e., aaa and ae are meas- 


ures of sodium and water excretion respectively, the origin corre- 
sponding to a condition of no reabsorption. 


A. Reduction of tubular Na’ reabsorption. 


It may be seen (Fig. 6) that for any given trajectory, reduced 
sodium reabsorption, due to decrease of Ena, favors increased 
water excretion. This effect is seen with the mercurial diuretics 
(Pitts and Sartorius, 1950; Duggan and Pitts, 1950; Smith, 1951; 
Farah and Koda, 1955), ‘Diamox’ Friedberg, Halpern, and Taylor, 
1952) and the xanthine derivatives (Sinclair-Smith, Kattus, Genest, 
and Newman, 1949). It is also seen that where Ewa is low, water 
reabsorption will be altered to a much less extent by changes of 
Cy than when Eng is high (cf. Grossman Weston, Borun and Leiter, 


1955). Similarly, if C,, is low, change of Ewa has comparatively 
little effect (Figs. 7 and 8). 


B. Reduction of tubular water reabsorption. 


Assuming that water diuresis causes a reduction of neurohy- 
pophyseal ADH activity (Verney, 1954) there will be a fall of C,,; 
Ena increasing. At the onset of a water diuresis the steady state 
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Na REABSORPTION —e— 


N. 


| where Es 20mV 


WATERs 
REABSORPTION 


where E=7OmV 


FIGURE 7, Indicating that where Eyq is low (= 20 mV) the change of 
water reabsorption, indicated by the black columns, caused by a given 
change of C,, is less than where Fg is high (= 70 mV). The arrows on 
the dotted lines represent decreases of Cy, i.e., movement from the left- 
hand trajectory to the right-hand one. 


Na REABSORPTION <= 


Oo 


f | where C,=O-OIC,, 


WATER 
REABSORPTION 


where C,=SOC,, 


FIGURE 8. (Water diuresis.) Indicating that where Cy is low (the 
right-hand trajectory), a change of Zwq has less effect upon water reab- 
sorption, indicated by the black columns, than where Cy is high (the left- 
hand trajectory). The arrows on the dotted lines indicate decreases of 


E'Na from 70 to 20 mV in each case. 


moves from A to B (Fig. 9) with an increase of water excretion. 
There can also be a tendency for sodium and chloride excretion to 
rise as A. R. Cushny observed (1917, p. 135). The sequence of 
events in hypotonic saline infusion is.similar (Cole, 1956). 
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C. The infusion of isotonic saline. 


When isotonic saline is infused there is a decrease of C, due 
to changes of body fluid volume although Ewa is not altered ini- 
tially 4 > B (Fig. 10). This will increase excretion of water and 
also_that of Na (Welt and Orloff, 1951; Leaf and Mamby, 1952). 
Reference to Figures 3 and 4 shows that these factors tend to 


Na REABSORPTION Sa 


WATER 
REABSORPTION 


FIGURE 9. The dotted line, AB, represents decreasing Cy with Eng 
constant. As the steady state moves from A to B, water and sodium re- 
absorption decrease, 


make the urine hypotonic with respect to plasma. As the urine is 
hypotonic and the infusion isotonic, the body fluids will tend to 
hypertonicity. Supposing this to cause a decrease of SRH secre- 
tion, Eva will decrease and the steady state will move from B to 
C (Fig. 10). Again from Figures 4a and 4b a reduction of Ena ata 
given value of C, will tend to cause a hypertonic urine in the 
steady state. Under these conditions, therefore, one would antici- 
pate the excretion of isotonic saline to commence with a phase of 
hypotonic urine excretion, succeeded by a period of increasing 
hypertonicity. The hypotonic phase has been recorded by Dastre 
and Loye (1889), Magnus (1900), and Green and Farah (1949). In 
1916 J. G. Priestly observed both a hypo- and hypertonic phase, 
as have Green and Farah (oc. cit.) and Cole (loc. cit.) (Fig. 11). 
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Na REABSORPTION 


WATER 
REABSORPTION 


FIGURE 10. (Isotonic saline diuresis.) The line AB represents a de- 
crease of Cy caused by reduction of ADH; BC represents a later phase 
during which #'y, decreases and C, remains unaltered. 


TIME hrs. 


Infusion of isotonic saline to 


Ris edt Se -Lo ml/min 


FIGURE 11. Isotonic saline diuresis. Ordinate represents urine so- 
dium concentration in mE g/L; abscissa represents time in hours. (Pooled 
results from 8 animals.) Note the initial phase during which the urine is 
hypotonic and the later phase of hypertonicity.- 
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IV. Conclusions. 


This theory implies that the action of the hormones is concerned 
with altering the steady-state reabsorption. From a physiological 
standpoint the shortcomings of the present treatment will be ap- 
parent: the hormonal control is certainly over-simplified and tubular 
morphology disregarded, Nevertheless, a number of physiological 
findings can be accounted for, and it is not unreasonable that the 
behavior of this model system may have some relevance to prob- 
lems of renal excretion of water and Na. 

It must also be emphasized that the treatment of the equations 
defining the model has been confined to steady-state solutions. A 
complete solution of these equations is being sought with the aid 
of an electronic analogue computer, and it is also hoped to extend 
this analysis to allow consideration of functional differentiation 
along the renal tubules. 
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OUTLINE OF A POSSIBLE MATHEMATICAL APPROACH TO 
THE PROBLEM OF THE EFFECTS OF ENVIRONMENT 
UPON THE INCIDENCE OF SOME PSYCHONEUROSES 
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In previous papers we have studied the distribution of the differences 
of two desires or inner drives in a group of individuals, with a view 
to developing a theory of mass behavior. The same concept is now 
applied to the problem of incidence of neuroses. - It is natural to con- 
sider that the incidence of neuroses in a group of individuals, which 
is characterized by a near equality of two opposite drives, will be partic- 
ularly high. Applying the results of a previous paper (Bull. Math. Bio- 
physics, 18, 323=36, 1956) in which the distribution of certain types of 
opposite drives has been connected with the distribution of incomes, ex- 
pressions are now derived which give the relative incidences of some 
types of neuroses for different income groups. The expressions derived 
hold only for certain types of assumptions and cannot, therefore, be di- 
rectly applied to actual data. They suggest, however, new possibilities 
in experimental and observational approaches to the problem. 


ie 


Landahl’s theory of conflicting stimuli (Landahl, 1938), which 
was originally developed in connection with a theory of psycho- 
physical discrimination, has found a very wide range of applica- 
tions (Rashevsky, 1948; 1951). Amongst other possibilities it 
offers a basis for a mathematico-biophysical approach to the prob- 
lem of some neuroses. 

Landahl studies basically the situation in which the excitation 
- of some center results in given reaction, while the excitation of 
another center results in a different reaction. If the two centers 
cross-inhibit each other, then the two reactions become mutually 
incompatible. » 

Such a model has the following property: If the intensity of ex- 
citation of one of the centers is sufficiently greater than that of 
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the other (situation A), then the reaction corresponding to the first 
center will be produced. If the two intensities, which may symbolize 
the intensities of two opposite drives, are nearly equal (situation 
B), then, depending on the neurophysiological constants of the sys- 
tem, either no reaction occurs at all, or the two reactions alternate 
at random with approximately the same frequency of incidence. 

Even in this highly oversimplified and unrealistic case of two 
simple reactions only, we may discern a certain remote analogy to 
some types of actual behavior. In situation A, the behavior may be 
said to be well integrated, with the organism ‘‘knowing what it 
wants.’’ In situation B, the organism’s behavior oscillates between 
two mutually contradictory activities, or the organism does not re- 
act at all. This may be considered as a very primitive case of 
neurotic behavior, and it occurs whenever two conflicting ‘‘drives”’ 
of approximately equal intensity are present. As the difference 
between the intensities of conflicting drives increases, the be- 
havior, though still oscillating between two contradictory activities, 
changes quantitatively. The frequency of the activity that corre- 
sponds to the stronger drive increases at the expense of the fre- 
quency of the other activity. While, for example, for equal con- 
flicting drives, both activities were equally frequent, now the one 
which corresponds to the stronger drive may occur in 75% of the 
cases, the other in 25%. If we ever had in reality such a simple 
situation as this, then the relative frequency of the two activities 
could be used as a measure of the neurotic state. A normal state 
would be characterized by one of the reactions occurring almost 
always; the completely neurotic state by a fifty-fifty occurrence 
of both. 

Actually, the situation is much more complex. Even the simplest 
activity of any organism is never a simple neurophysiological reac- 
tion of areflex type. It consists of a very large number of coordinated 
individual reactions. When an animal exhibits two contradictory 
activities, for example, taking food and running away from it in 
fright, these two activities may still have some elements in com- 
mon. Thus, the use of the extremities in locomotion either toward 
the food or away from it involves essentially the same elementary 
neuromuscular reactions. In an actual situation, therefore, an ani- 
mal may, under conditions of neurosis, exhibit neither one nor the 
other of the two activities, but only that constellation of reactions 
which is common to both activities, but which in itself does not 
constitute an organized pattern. 
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This complicates the theoretical problem considerably, but leaves 
its essential features unchanged. Whereas a single drive results in 
a coordinated reaction pattern, two conflicting drives of nearly equal 
intensity result, in general, in an uncoordinated set of reactions. 

In Landahl’s original approach to the problem (1938), the model 
was applied to psychophysical discrimination. 

If, for example, an individual is given two weights, one 400 
grams and the other 200 grams, and is asked which is the heavier, 
he will without mistake give the correct qgnswer. If, however, the 
weights are 400 and 405 grams, in almost 50% of his trials he will 
give the wrong answer or will not be able to give an answer at all. 
Landahl’s equations permit us to calculate, for any given pair of 
weights, the percentage of correct answers. The calculated values 
agree amazingly well with the observed ones. 

Such an experiment may have nothing to do with neuroses. Even 
if we are off 50% in our judgments we do not necessarily exhibit 
neurotic behavior. If, however, the subject would be told that his 
life depended on always making a correct judgment, we probably 
will agree that he may develop neurotic behavior almost before the 
experiment starts! It was a similar experiment on dogs that led to 
the discovery of conditioned neuroses in animals. 

The next task of the theory would be to explain the process of 
‘*conditioning’’ or production of neurosis. What is the neural mech- 
anism through which certain environmental factors result in the 
formation of equal conflicting drives? This leads us to the mathe- 
matical theory of learning. Again, some quite successful results 
were obtained in the theory of some types of learning. We say 
‘some types of learning,’’ because it seems to be highly unlikely 
that all learning is mediated by the same mechanism. Learning to 
type and learning to memorize the Iliad are both in the ultimate 
analysis reducible to the same laws of neurophysiology, and the 
neural mechanisms involved are built of similar neurons, Never- 
theless, the two mechanisms are probably as different as a tele- 
phone exchange and an electric-power plant, though the operation 
of both is ultimately reducible to the same laws of electrodynamics, 
and both use wires, coils, magnets, etc. The mathematical equa- 
tions derived for the rate of learning of some tasks and verified ex- 
perimentally are likely, with some modifications, to be applicable 
to the rate of formation of some neuroses in experimental situations. 

By considering the conflicting stimuli to be of endogenous na- 
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ture, we arrive at a theory of behavior in situations in which an in- 
dividual is confronted with two inner desires, or drives, towards 
two different and mutually exclusive behaviors (Rashevsky, 1951). 
Here again, by using the same argument as above, we may expect 
an individual to behave neurotically if his two opposite drives are 
equal, or nearly equal. This is also in line with psychiatric findings. 

Other types of conflicting stimulations, of a more complex nature, 
are known and are also conceivable neurobiophysically. One of the 
two conflicting stimuli may be endogenous, of the nature of an in- 
ner drive. The other may be exogenous. A situation of that nature 
may occur, for example, when an individual craves to perform an 
activity, but when external circumstances make him perform an op- 
posite activity. Such ‘‘frustration’’ situations should theoretically 
also lead to the development of neuroses and, as is known, they 
actually do. 

Thus any environment which offers contradictory situations is 
more likely to lead to the development of neuroses than an environ- 
ment which does not offer such situations. 


Il. 


One important factor which determines to a large extent the en- 
vironment of an individual is his economic status. Whenever the 
economic needs of an individual exceed his actual income, situa- 
tions of a conflicting nature are likely to occur. If z denotes the 
income of an individual and 7 his need, in monetary terms, then 
the greater z’.- z, the more likely and more frequent will be various 
conflicting and frustrating situations. Putting 

visa, (1) 
we may Say that, in a group of individuals characterized by a given 
~ and a given 2, the frequency of occurrence of conflicting situa- 
tions, and therefore also the frequency of occurrence of neuroses, 
will be a function f(2,z), such that df/dz > 0. As regards the char- 
acter of dependence of f on 7, we have no direct indications about 
it. It would perhaps seem natural to expect that f(i,x) will in- 
crease with increasing ratio z/i. Psychologically this would mean 
that an individual, who can satisfy more of his desires than can 
another, will be less likely to feel frustrated by the impossibility 
of satisfying a given number of desires than will be the other in- 
dividual, However, much more complicated relations are possible, 
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and perhaps even likely. In a first, abstract, theoretical approach, 
in which more emphasis is laid on the methodology rather than on 
concrete results, we are justified in choosing the simplest possible 
assumption about f(z,z). Such an assumption, which may be a crude 
first approximation, is to make f independent of i, and a function 
f(x) of 2 only. This is what we shall do in this paper as a pre- 
liminary step to a subsequent study of more complex and more real- 
istic situations. 

As in a recent paper (Rashevsky, 1956), we shall consider a 
distribution function which gives the fraction of individuals in a 
society that are characterized by an income between i and i+ di 
and by a need between 2’. and 7’-+ di”, Denote this fraction by 
G(z,2’) didi’. We have 


[ I G(i,i’) didi’ = 1, (2) 
ae 


where z,, and 7 are respectively the smallest occurring income 
and the smallest need. If G(i,2’) approaches zero sufficiently 
rapidly, when z — 0 and z’.— 0, then without appreciable error we 
may use zeros as lower limits in (2), 
The expression, 
Ni(é’)di’-= ai GB(i,v’) di, (3) 
tm 
represents the distribution function of needs in the population. 
The expression, 


Na(2)di = di [ G(i,i’) di’, (4) 
vn 
represents the income distribution. The distribution No(Z) is readily 
measured and known for many populations. The distribution N ,(2’) 
could be and in some instances has been measured, at least roughly, 
by appropriate questionnaires. By a questionnaire method the func- 
tion G(z,i’) can also be measured, » 

If No is the total number of individuals in the population, then 
the number of individuals whose income lies between 7 and 7 + di 
and whose needs exceed their income by an amount between « and 
2+ dz is given by : 


NoG(i,i + 2)didz. (5) 
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Of these individuals, the fraction f(z,7) will develop neuroses. 
Hence the total number of individuals whose income is between z 
and i + di and who develop neuroses is 


n(i)di = Nodi ic Gli,é + 2) f(é,2)de. (6) 
0 


Therefore, the fraction v(z) of neurotics amongst the individuals 
with income between i and 7 + d is equal to n(z)di divided by the 
number of individuals who have an income between z and z+ d. 
But that number is equal to NoN2(z)di, where Noa(z) is given by (4). 
Hence, from (4) and (6): 


G(i,i + x) f(2,2)de 
v(z) = 2 ______— ; (7) 
iE Git”) ai’. 


m 


Thus, from the knowledge of the functions G(z,2’) and f(z,z), we 
obtain the distribution of incidence of neuroses as a function of 
income. This distribution of incidence v(z) is expressed in terms 
of the parameters of G(z,z’) and f(z,z) and in terms of 7},. 

The fraction v; of neurotics in the whole population is given by 


co 0 «8 


vr =| i v (a)G(i,2’) didi’. (8) 


and thus is also expressed in terms of the parameters of G(z,i’) and 
f(i,z), as well as in terms of 7, and 7%. 


Ill 


The whole problem has, however, another entirely different as- 
pect. In a previous paper (1956) we discussed the effect of the 
distribution G(z,z’)} upon the distribution of the tendencies of the 
individuals toward the following two opposite and mutually exclu- 
sive behaviors: the one is the preservation of the existing socio- 
economic structure and the other, its change, regardless in which 
direction. The argument was as follows. 

An individual, for whom ¢ > 2’, is likely to be satisfied with the 
existing situation and tend to support its preservation. This tend- 
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ency will be the greater, the greater i ~ i’, On the contrary, an in- 
dividual with 7 <2’ will have psychological grounds to be dissatis- 
fied with the existing situation and try to change it. The desire 
for such a change will be the greater, the greater 7” — 7, Indi- 
viduals with z = ¢” would be just satisfied and are likely to be in- 
different one way or another about any changes in the existing 
situation. Therefore, 


b=i-i. (9) 


' may be considered as a measure of the tendency of an individual 
towards either the status quo (¢ > 0), or towards a change (¢ < 0). 
When ¢ = 0, the individual is indifferent. The quantity 4 thus rep- 
resents a measure of a tendency towards one of two conflicting 
behaviors. 

Let N(¢)ddé denote the fraction of individuals whose ¢ is be- 
tween ¢ and ¢ + d¢. As has been shown previously (Rashevsky, 
1956) N(¢d) can be expressed in terms of G(z,z’) in the following 
manner: 


For ¢ > 0: N(#) = N;(¢) = [ Qli-+ d, 0) di"; (10) 
0 

For ¢ <0: N(¢) = N_() = ‘i G(i,i — )di. (11) 
(0) 


In deriving equations (10) and (11), the quantities 7, and tn 
were considered as zero, which does not introduce any restrictions 
in principle. We shall throughout this paper also put tm = 7, =0 
for simplicity. 

It is readily seen that, although composed of two different branches, 
N,(¢) and N_(¢), the function M(¢) is continuous at 6=0. We 
have N,(0) =N_(0). The first derivative dN(¢)/d¢ is also con- 
tinuous at ¢ = 0. 

According to the picture outlined in the introduction to this paper, 
we would expect that individuals with a sufficiently large positive or 
negative ¢ will exhibit an integrated normal behavior. But individ- 
uals with ¢ = 0, or with very small positive or negative ¢, will be 
exhibiting an erratic behavior characteristic of neuroses. The nar- 
rower the interval A around zero, in which ¢ lies, the more random 
or erratic the behavior of an individual will be. A completely er- 
ratic behavior will be exhibited by individuals with d= 0. In any 
society, however, it is a matter of more or less arbitrary conven- 
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tion or tradition, or medical custom, as to how much the behavior 
of an individual must deviate from ‘‘normal,’’ in order to consider 
the individual as being ‘‘abnormal.’’ Hence the quantity A is to 
the same degree arbitrary and reflects the relativity of the notion 
‘abnormal behavior.’? Thus abnormal behavior will be exhibited 
by individuals with 


-tA<<+ HA. (12) 


It may seem at first rather absurd to assume that neurotic be- 
havior will be exhibited by individuals whose incomes are equal to 
their needs. One would a priori expect perfect contentment in such 
individuals. The seeming absurdity of the assumption, however, is 
at least considerably reduced when we keep in mind not merely the 
role of conflicting stimuli in producing neurotic behavior but also 
the nature of the stimuli and of their responses. The two conflict- 
ing stimuli may be, for example, the desire to swim and the desire 
to walk. No one is likely to become neurotic in this situation. 
When the two desires are equal, the individual will simply on the 
average spend half of his available time in swimming, half in walk- 
ing. No one would call such a behavior erratic, and no particular 
mental stress is attached to it. Any time the individual is tired of 
swimming he may, if he chooses, take a walk, and vice versa. Let, 
however, an individual have equal desires to follow two mutually 
exclusive carriers, each of which will take up all of his time. The 
result may be either a prolonged indecision, which is quite likely 
to result in mental stress and neurotic symptoms, or the individual 
will start on one carrier, abandon it after a while for the other, 
then again return to the first, etc. His behavior will justifiably be 
called erratic in this case, and his achievement will be nil. Neu- 
roses on such a basis are well-known. » 

Even a conflict between swimming and walking, or rather run- 
ning, may lead to neuroses, in appropriate conditions. Imagine the 
following situation which may appear somewhat farfetched but is 
quite possible. A man is attacked by a bear, which is equally dan- 
gerous on land and in the water (Sanderson, 1956). The man may 
either run, or jump into the water and swim. If the indecision does 
not cost him his life, it certainly is likely to result in severe 
neurosis, 

The quantity 4, as defined by equation (9) is a measure of the 
desire toward the following specific set of two mutually exclusive 
activities: to preserve the economic status quo of the society or to 
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change it in any manner. Most individuals with ¢ =0 will be merely 
indifferent to the whole problem. Also by far not all of the individ- 
uals with ¢ =2 — 7’-<0 will have a desire to act towards an eco- 
nomic revolution. Nor will all individuals with d > 0 constantly 
worry about the preservation of the status quo, But among all 
those individuals there will be a fraction, possibly but not neces- 
sarily a small one, which may be called ‘‘socially conscious.’’ If 
such a man has a ¢ < 0 he is likely to take an active part in reform 
and revolutionary activities. If he has a ¢>0, he will actively 
oppose such activities. If he has a ¢ = 0, he may well be torn be- 
tween two opposite drives to the extent of developing a neurosis. 

Therefore not all of the individuals with a ¢ satisfying inequality 
(12) will necessarily develop neuroses, but only a fraction a of 
them, which are strongly socially conscious. 

The type of behavior with respect to which conflicting situations 
arise according to the mechanism discussed in section II of this 
paper is more likely to be of the nature of the individuals’ daily 
routines. The behaviors which are involved in the mechanism dis- 
cussed in this section are in the field of large scale social activ- 
ities. There is nothing surprising, therefore, that neurosis may be 
produced in the two cases by opposite factors. In the first case a 
large z’-— z is favorable to the production of neuroses, in the sec- 
ond a value of z’--- 7 near zero is essential. Both factors are likely 
to contribute their share to the total incidence of neuroses. And 
there can be no doubt that there are numerous other ways in which 
the economic status can affect mental health. Some of those fac- 
tors are to be considered in subsequent publications. In this first 
orienting study we choose a couple of the simplest abstract situa- 
tions for a mathematical approach. Only through such a study first 
of the simpler, then of the more complex factors may we hope to ar- 
rive in the future at a more comprehensive theory which will enable 
us to analyse and interpret actual facts. 

In general, the fraction % of the socially conscious individuals 
may be a function o(¢) of d. Putting (0) = 40, we thus see that 
the total fraction vj of individuals in a society, characterized by 
the distribution function N(d) [equations (10) and (11)], which ex- 
hibit neurotic behavior is given by 

+74 
ier ee i N(¢)dd. ; (13) 
1 


ZtA 
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When A is sufficiently small, as is likely to be the case, we may 
with sufficient approximation write, instead of (13): 


v7 = 40 AN(0). (14) 


We may now ask for the incidence v’{z) of neurotics amongst in- 
dividuals whose income is equal toz. The total number of neurotics 
with income between z and 7 + dz is 

i+FA 

AoNodi G(i,0’) ae’; (15) 

ita 
for this is the number of individuals whose income lies between 7 
and i+ di and whose need is between i-}A andi+ 3A. The 
total number of individuals with an income between z and z + d is 
given by No times the expression (4). Again putting 7, =0 and 
dividing the number of neurotics given by (15) by the number of in- 
dividuals NoNo(2)di, we find: 
i+tA 

G(t,2’) ar’ 

* 1 
; i-FA 
v°(z) = Xo : (16) 


co 


f G(i,é) ai? 
0 


If A is sufficiently small, we may use the same approximation 
which leads to (14). Then instead of (16) we obtain: 


ao A G(i,i) 


i Gi,” ai? 
0 


If the factors discussed in this section and in section II were 
the only operating ones, then the incidence of neuroses, due to 
both, would be given by vr +v‘, respectively by v(%) + r“(4), 
where v7, vr, v(z), and v’(2) are given by (8), (13), (7), and (17). 

The situation discussed in this section holds for the case of ab- 
sence of immitative behavior (Rashevsky, 1951, 1956). The effect 
of the latter, in its simplest form (Rashevsky, 1956), is to shift 
the curve N(¢) to the right or to the left by an amount w*, deter- 
mined by equations discussed previously (Rashevsky, 1951, 1956). 
This results in a ‘‘uniformization’’ of the behavior of the whole 


v’(z) = (17) 
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society, a large majority exhibiting now one of the two mutually 
exclusive behaviors, and only a small minority—the other. Inas- 
much as for physical reasons N(~-) = N(+0) = 0, such a shift will, 
if sufficiently strong, result in a decrease of v7, as given by either 
(18) or (14), because M(+yW*) is small for sufficiently large *. 
The actual calculation of v7 proceeds as follows: Suppose the 
shift is in the direction of positive ¢’s, that is, for a ‘‘conserva- 
tive’’ behavior, and let it be given that a fraction x <q<1 of all 
the individuals exhibit the conservative behavior. Let P1(¢) de- 
note the frequency of conservative behavior exhibited by an indi- 
vidual characterized by a given value of 4, (Rashevsky, 1951, p. 
88). Then the difference between the number X of individuals ex- 
hibiting ‘‘conservative’’ behavior and the number of individuals, 
Y, exhibiting the opposite behavior is given by 


+00 


beens a Bf N(¢) [2P1(¢ + ¥*) - 1] d¢. (18) 
We also have X + Y = No. Hence the fraction g is given by 
* en | a 
g= ae Ae =f N(¢) [2P (¢ + $*) - 1] dd. (19) 
) een 


The right side of (19) is a function of y*. Equation (19) may be 
solved for w*, giving the latter as a function w*(g). If the NM(¢) 
curve has been shifted by the amount y* to the right, then the 
value of its ordinate at ¢ = 0 is now M(—y*). Introduced into (14) 
this gives 


vp = do AM(-Y*). (20) 


The more exact expression which corresponds to (13) is obtained 
in a Similar manner. 


IV. 


Just to illustrate the type of relations which may be obtained by 
an explicit evaluation of the expressions derived in the two pre- 
ceding sections, we shall choose as an example a distribution func- 
tion G(i,2’) which has been considered in a previous paper (Rashevsky, 
1956). 

If ¢ and 2’ are not correlated, then G(i,z’) becomes simply the 
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product of the distribution function N;(2’) of the needs and N2(z) of 
the income. Thus 


G(i,t’)didi’-= N 1(2’)No(t)dr de. (21) 
For N,(2’) we shall choose, for illustration purposes only, (Rashev- 
sky, 1956): 


bP - 
NCE et? eae (22) 


where 6 is a constant. For No(z) we choose 


No(2) = ce~*?, (23) 


where c is a constant, 
¥ quations (10) and (11) which define N(¢) now become (Rashevsky, 
1956, p. 328): 


bee 
For ¢ > 0: NBT SING) ae io (24) 
63 1 +2 
For 6 <0: Nd) = N_(¢) = > ase = os + aa eoP, (25) 


For f(z,z) we make the simplest possible assumption, namely, 
that it does not depend on 7 and is proportional to z. in other words, 
denoting a constant by a, we put 


fl(e,2) = az. (26) 


Expression (7) now becomes: 


oo 


oleeas [ aNi(é + 2)de. (27) 
0 
Introducing (22) into (27) we obtain: 
;, Dirt | hv and 
v(i)=@ ee + 2i4+ ;) ie (28) 


By differentiating (28), we find that dv(i)/di <0 everywhere. 
Thus v(z) decreases monotonically from its greatest value 3/6 at 
¢=0. In other words, with the specific assumptions about N;(7’), 


No(z), and f(z,z), the incidence of neuroses will be highest among 
the poorest people. 


ENVIRONMENT AND PSYCHONEUROSES 53 


Expression (8) gives (with 7%, = im = 0): 


co 


a. a v (4)Na(i)di. (29) 


Introducing (23) and (28) into (29) we find: 


6b2c + 8bc? + 3c3 


noes el OS aseyey |, 2° ) 


which expresses v7 in terms of the two parameters 6 and c of the 
distributions N;(2’) and No(z). The derivative of vr with respect 
to c is equal to 
dvr 6b3 + 462¢ + bc? 
= a ee ee eee 
de 6(5 + c)* 


>, (31) 


Thus vr increases monotonically with c. For c = 0, vr = 0; for 
3a sas 
C=e,yr= axe But 6 and c are related to the average need 7 and 


the average income z by the expressions: 
3 ae 
= 3B? = oa (32) 


if the distribution functions are given by (22) and (23). Hence we 
see that vy increases as the average income 7 decreases; and as 7 
tends to zero, vr tends to ai’. This, of course, holds only for the 
particular choice of the distributions N1(2’) and No(z) made above. 

For the derivative of v7 with respect to 6 we have: 

d 125% ¢ + 18b7.c? + 125c* + 3c* 
at. -g ——_____—__,__— <0 (33) 
db 67(b + c)* 

Thus vr decreases monotonically with 6. For } =0, vr =»; for 
b=00,vp=0. Or in terms of 7: vr = for 7” =~; and vr = 0 for 
7’= 0, 

We shall now give the expression for v7, using the approxima- 
tion (14). From both (24) and (25) we see that 


3 


(hl) ere 34 
Hence, introducing (34) into (14), we find: 
bo A b%c 
Vege (35) 


th.+ ae 
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For c= 0 and c=~, vy = 0. Hence it has at least one maximum. 
It is readily seen that only one maximum exists. We have 


ee hg Nhe (36) 


nax = = SS at 
pee G7) 
Introducing (37) into (35) we find the maximum value of v7: 
4bao A 440A 
= =o, 38 
"Tmax ~ 97 97” (22) 
We also have 
dv’ 340 A b2c? 
ws : > 0. (39) 


‘db (b+0)4 


This shows that vr increases monotonically with 6, from vr = 0 
for 6 =0, to v7 = do Ac for b=. In other words, »’p decreases 
with 7’. 

We now derive an expression for v’(z). Introducing (21) into (17) 
we find: 


v’(i) = do ANi(A), (40) 
or, because of (22): 
Xo Ab ; 
y(t} =<— ao (41) 
The value of v’(z) has a maximum for 
2 9% 
Chak. =o eee e 


and its value at the maximum is: 


2aoAb 6aAoA 


a ee (43) 


v(t) max Se 
e 
where e is the basis of the natural logarithms. 

Equation (38) expresses the maximum possible value of the total 
incidence of neuroses in a given society which is characterized by 
the distributions N ;(2’) and Na(z), given by (22) and (23), and which 
has a given average per capita need, 7’. Equation (41) gives the 
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distribution of the incidence of neuroses for different income classes 
in the same society, and equation (43) gives the incidence for the 
income class in which it is highest. 


V. 


We shall now discuss the possibilities of empirical verification 
of some of the above conclusions, as well as how they can be 
used for the interpretation of some observational data. 

First of all we may point out one important practical conclusion 
which must be kept in mind in all the statistical studies of the in- 
cidence of neuroses in different income groups. This incidence 
depends not only on the income but on the distribution of the in- 
comes and of the needs in the whole society. This conclusion does 
not depend on the particular choice of the distribution functions 
made in section IV of this paper. It follows directly from equations 
(7) and (17). We cannot, therefore, compare the incidences of neu- 
roses in two different income groups, if the two groups are taken 
from different communities in which the distributions of incomes 
and needs may be different. 

If z and z’ are not correlated, a situation which is unlikely, then 
G(i,z’) is given by (21). The income distribution N2(z) can be de- 
termined and has been determined for many contemporary societies, 
The distribution N ;(z’) of needs may be more difficult to determine 
because there is no objective yardstick for measuring it. How- 
ever, by a proper questionnaire method it is possible to gain an 
idea of N;(z’), expressing z’-in monetary units. By a combination 
of direct economic census with the use of appropriate question- 
naires, it should be possible to determine G(z,z’) in the general and 
more likely case where 7 and 2’ are correlated. 

Once we have G(z,i’) and express it by an adequate empirical 
equation, we may make different plausible assumptions about 
f(i,z) and for each assumption derive the values of v, and the 
function v(z), as well as vy and (i). Those then can be checked 
again by statistical observations. If the two factors, which con- 
. tribute to the incidence of neuroses and which are discussed re- 
spectively in sections III and IV, are the only present (an assump- 
tion which again is not likely), then it is vr + v7 and v(z) + v’ (2) 
which would be given by statistical observations. The assumption 
about f(z,x), which leads to expressions for vr + v7 and v(?) + v’(z) 
that agree best with observations, would be considered as the one 
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corresponding best to the real situation. Once such a form of (2,2) 
is found, it may, in principle at least, be directly compared with 
observations. If we had available all the desirable statistics on 
mental patients, then it would be possible to determine directly 
what percentage of persons with a given income 7 and given need 
v’.=i+ x suffer from neuroses. 

A different approach would be to derive G(z,z’) from some other 
theoretical considerations, as a function with arbitrary parameters. 
Then, using the methods illustrated in the examples of section IV, 
we could derive such expressions aS v7,,,,5 OF v(¢)max, OT max, all 
in terms of the arbitrary parameters. By determining vrs VCe)maxs 
and tmax from statistical observations, we may then calculate the 
values of the parameters. Those calculated values could then be 
checked against observations by a direct determination of G(z,2’). 

All the above program is actually much more easily suggested 
than executed. But this is the case with every exact science. The 
important thing is that the different conclusions of the theoretical 
approach outlined here can be verified empirically in principle. 
Therefore they do suggest different empirical approaches which 
otherwise might not have been thought of. The possible verifica- 
tions are, as we see from the above, of the following type: we de- 
rive the values of some constants, or the shape of some function, 
from a set of data of different type and see whether the findings 
thus obtained agree with each other. 

A more detailed psychiatric analysis of the patients may eventu- 
ally permit the separation of the empirically found sums v7 + uP 
and v(z) + v’(z) into their components. As we have seen in section 
III, (p. 49), the nature and possible etiology of the neuroses, the 
incidence of which is expressed by v7 and v(z), is somewhat dif- 
ferent from the nature of the neuroses, the incidence of which is 
expressed by vy ‘and v’(7). 

Or again, we may use the relations discussed in this paper to 
determine N ;(z’) from a known No(Z) and f(z,2). 

As said above, the two factors producing neuroses discussed in 
sections II and III are not the only possible ones. Considering 
other factors, we shall complicate the situation further both mathe- 
matically and empirically. But as has been the case in older ex- 
act sciences when we have a complex situation with a large num- 
ber of factors involved, mathematical analysis of the abstract, 
oversimplified cases that involve each only one or a couple of 
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factors helps us in the long run to separate those factors empiri- 
cally by pointing out the way to such measurements which can be 
used to effect such a separation. 

The following very important point must be emphasized. While 
the economic status is undoubtedly an important factor in deter- 
mining the incidence of some mental diseases, in particular of 
Some neuroses, it is never the only factor. Of two individuals 
with equal economic status and equal needs, one may develop a 
neurosis, the other not. The first may have reached his present 
status through struggling from a lower one in which conditions 
favoring the production of a neurosis were abundant. The second 
may have been free of economic worries and conflicts all his life. 
Purely non-economic factors also enter the picture. 

To take into account mathematically this situation, we must com- 
plicate the problem still further by considering effects of social 
mobility upon the incidence of mental disease. The problem be- 
comes much more complicated not only mathematically but also 
empirically. A mere statistic of incidence of mental disease in 
different economic groups now becomes meaningless. The diffi- 
culties of experimental approach are augmented many times. But 
this again is the way of all science, and again it is in just such 
cases where an intentionally oversimplified mathematical analysis 
is likely to help us to separate the individual factors and to unravel 
their mutual interactions. As someone once said, nature is under 
no obligation to scientists to be simple and easily understandable. 
If we are afraid of present or future complexities, we may just as 
well abandon any attempt at a scientific approach and progress. 

By way of an outline of further theoretical studies in the direc- 
tions indicated here, we may mention the following problems: 

First of all a more detailed theoretical study of the function 
G(i,i’) should be made, thus narrowing the possible range of its 
shape. Different mechanisms of correlation between 7 and 7’. are 
conceivable, and they should be studied mathematically, beginning 
with the simplest ones. 

The theory outlined in section III should be amplified by intro- 
ducing finite thresholds for the effectiveness of conflicting stimuli, 
as has been done by H. D. Landahl (1938, 1950). The analogy be- 
tween the behavior under the influence of conflicting stimuli for small 
values of 4 becomes then much more pronounced. (See section I.) 

The considerations of a previous paper (Rashevsky, 1956), which 
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lead to the expression for N(¢) in terms of G(z,z’) Lequations (10) 
and (11)] may be applied, mutatis mutandis, to distributions of 
social status. An approach to the theory of such distributions has 
been made previously (Rashevsky, 1951). The theory indicates a 
connection between the distribution curves and various other socio- 
logical parameters, such as the frequency of contact between in- 
dividuals of different social standing. Thus a number of other socio- 
logical parameters will be introduced into the expressions for in- 
cidence of neuroses, and in this way a number of other suggestions 
for experimental approach will be found. 

A complete parallelism between theoretical and experimental ap- 
proach is the ideal of any exact science. This ideal, however, is 
rarely met in reality. Frequently, experimental research precedes 
by too much the theoretical approach. Valuable results may be 
found experimentally or observationally which are, however, not 
too well adapted to an interpretation by a theory because no theory 
can be developed at once to fit a given complex experimental situ- 
ation. A mathematical theory has to start with a study of over- 
simplified, abstract cases, before it can reach a stage directly ap- 
plicable to a given complex actual situation. This may be con- 
sidered as a sort of generalization of the saying that there are no 
royal roads in mathematics. 

On the other hand, sometimes a theory is developed in consider- 
able detail, and it makes a number of predictions which are verifi- 
able. But it may take years until the experimental verification is 
forthcoming. Such was the destiny of some of J. C. Maxwell’s 
work, such was the destiny of the mathematical prediction of the 
law of transformation of matter into energy. 

In those sciences in which mathematical approach is still new, 
such situations are particularly apt to occur. Only patience can 
help overcome the difficulties thus created. 

To give an example which applies to the problem treated here, 
we may-mention that very valuable data on the incidence of mental 
disease in different socioeconomic classes have been obtained by 
A. B. Hollingshead and his co-workers (Hollingsnead and Redlich, 
1954, 1955; Redlich, Hollingshead, and Bellis, 1955; Hollingshead 
and Freedman, 1955). Though the data obtained by them are not 
given in the form of curves, they indicate definite quantitative re- 
lations, Naturally, a practical sociologist and psychiatrist may 
expect a theoretician engaged in applied mathematical studies to 
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interpret mathematically those empirical results. From the point of 
view of such an interpretation, the contribution of the present paper 
appears to be zero. However, the conditions of the experimental 
study, which may have been the simplest possible one for an em- 
pirical research, do not happen to be the simplest ones for a theo- 
retical approach, and vice versa. It is our purpose in future theo- 
retical studies to approach the conditions which were characteristic 
of the empirical studies. A cooperation on the part of the empiri- 
cal researcher, which perhaps may require some modification of 
the empirical approach, would be of great help to the applied math- 
ematician. 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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The methods of C. W. Sheppard and A. S. Householder (Jour, App. 
Physics, 22, 510-20, 1951), H. D. Landahl (Bull. Math. Biophysics, 16, 
151-54, 1954), and H. E. Hart (Bull. Math. Biophysics, 17, 87-94, 1955) 
are extended to non-conservative, non-steady-state systems whose com- 
partment volumes are changing. It is shown that complete determination 
of the transfer, localization rates, and volumes is possible for systems 
in which there is access to all compartments. Information obtainable 
from experiments wherein some compartments are not accessible is noted 
for mammillary and catenary systems. If intracompartmental mixing is 
slow it is shown that steady-state systems can be interpreted in terms 
of non-steady-state systems whose transfer and localization character- 
istics are then functions of the time after and method of tracer adminis- 
tration. The importance of such intracompartment mixing is discussed. 
Types of steady and non-steady-state systems are outlined and the ap- 
plicability of various methods of analysis indicated. 


Introduction. Investigations of multicompartment systems have 
been concerned principally with steady-state conditions (Gellhorn, 
Merrell, and Rankin, 1944; Flexner, Cowie, and Vosberg, 1948; 
Solomon, 1949; Sheppard and Householder, 1951; Hart, 1955— 
hereinafter called reference I). It has been shown, however, that 
tracers can also be used to investigate non-steady-state, non- 
conservative systems (Landahl, 1954b). Since growth and pathology 
often involve variable rates and volumes, it may be of interest to 
outline in some detail tracer methods whereby such systems can 
be studied. It will therefore be a purpose of this work to consider 


tracer experiments in N compartment systems of which M compart- 


Research supported by the Atomic Energy Commission, Contract AT 
(30-1)-1551. ; 


61 


62 H. E. HART 


ments are accessible, where experimental errors are small and for 
which the volumes and transfer rates are continuous, ‘‘well-behaved’’ 
functions of time. 

In Section I, the complete solution? to such a system is obtained 
assuming access to all compartments and the application of the 
method to special cases outlined. In section II, the relation of 
intracompartment diffusion processes to non-steady-state systems 
is explored. In conclusion the types of steady and non-steady-state 
systems are briefly described and the applicability of different 
methods of analysis indicated. 


Table of Notation® 


z,,(¢) = concentration of freely diffusible labeled material of the 
Ith species in the zth compartment. 

V,(¢) = volume of the 7th compartment. 

r, jO = rate of fluid transfer from the jth to the 7th compartment and 
r,; 1s not defined. 4 

r,-(t) = total rate of fluid transfer from external sources to the 7th 
compartment. 

Te1%e1 = tate of entry of the /th species into the 7th compartment 
from external sources. 

r,,(t) = rate of excretion from and or localization in the zth 
compartment. 

—a,(t) = 7;,(¢)/V,() = rate of fluid transfer from the jth to the ¢th 
compartment per unit volume in the 7th compartment. 


N,e 
a,,(t) = D2 ry O/VO) where >.” indicates j #2. It represents the 
j i 

rate of transfer per unit volume in the ith compartment of fluids 
from all other sites. 

@(7,p) = average concentration of freely diffusible labeled material 
in compartment 7 of a steady-state system. 

a, ;(t,p) = effective transfer coefficient in a steady-state system for 
which intracompartment mixing is slow. + represents the time 
after injection and p the site of injection. 


2 Fy : 
Determination of rates and volumes from specific activity measurements. 
For a discussion of further interpretations of 2;,%;j, and V; see foot- 
notes on Page 88 of reference I. 
“Note that the subscri i i 
at the subscripts are interchanged relative to reference [3 in 
the present notation t,; takes the place of Tie 
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I. Continuous Measurements in All Compartments. 


If the volume V,(¢) and transfer rates ti j are assumed to be 
continuous functions of absolute time only$ (i.e., not dependent 


upon the time or site of injection), it follows from the definitions 
in the Table of Notation that: 


Gee ieedy 
Viti, + @,, V, = Sai (7, jen T%ia)> ae ee (1) 
j=1 jwileeN 
and 
N,e 
V, = 15 (4 =F) (2) 
j 
Let 
r= —ViG;; RS Pr ee Be ie ee See (3) 
N,e > 
@,,=- mee a,, where oe indicates j #7. (4) 
j j 
N,e 
Then @i.+ ac 0,25, = 0, (5) 
j 
N 
or &,, + a A; ,%;1 = ts 01% e1/V;: (6) 
j 
It follows that a;;= a; 25 a, 7/V, ’ (7) 
where 


Tiga My > Uy jpaiio NK 


Bit Mae Se SAS Ics (8) 


in°°* 2j—1,.n ~ Yin ®j41,N°°* *NN 


aR SoG = 1 Piet? el oy fo ome ee” 


fae oe ea, le (9) 


Ty 1 eee @3—1,NieN%eN~j+ LN Lae) InN 


5This is not true for diffusion problems; see Section II. 


64 H. E. HART 


and 7,.,a,, is the rate of entry of the Ith species into the zth com- 
partment from external sources. 
From the conservation of tracer material 


N N t 
ny V; oe ~> [ (7; 21% ¢1 a 7 oi% i) ats bw wedNz (10) 
b i 0 


L 


Since the system cannot discriminate between isotopes, r,, does 
not depend upon J. Differentiating equations (10), noting that 


N N 
"ei >. (r5- ya) + Tie — Vi fo ~~ Tiel? (11) 
j I 
and using equations (8) and (9), it follows that 


N N 
° + * + * 
S eat) (—V,a,; - a,5 + Via, + Qj )® x 
u j 
N N 


= Fa Ps Ca np Tem “il ’ (12) 


i m= 1 


where mis simply a dummy index. 

Equation (12) is a set of N linear equations in the N unknown 
functions V,. Since r,,, and z,, will probably be known, it follows 
that in principle equations (3), (7), (8), (11), and (12) can be used 
to obtain the general solution for a non-steady-state, N compartment 
system in which there is continuous access to all compartments 
and for which intracompartment concentration gradients are van- 
ishingly small. 

If, as will often be the case, all labeled material is injected 
simultaneously and instantaneously r;.,2,,=0 for ¢¥ ¢, and equa- 
tions (6) and (7) reduce to 


N 
& + i, G;,2 5) = 0. (138) 
j 
and 
epi aeetas E Ss cadet. 
eat eld : - 
of eke ee a a ead OF. emer | 


IN°** Toy ny ~ Vin Pai ness Cyy 
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Except for a change in notation, equations (14) are the same as 
equations (2) in the work of Landahl (1954b). If isotope 1 is in- 
jected in compartment 1, isotope 2 in compartment 2, and so forth,°® 
the diagonal terms (i.e., 2,;) in [exa are dominant initially and 
|X|=0. Since the denominator does not vanish in the early 
phases, it follows that a,;(t) will be uniquely determined in the 
beginning. Since a sufficient (though not necessary) condition for 
|X | = 0 is that the th column vector z,, be a constant multiple of 
the jth column vector © sys it seems reasonable that the a; (t) will 
be uniquely determined as long as in all compartments the relative 
distribution of the different injected isotopes continues to reflect 
their initial conditions or sites of entry. 

Although equations (8) and (9) are identical in appearance to 
equations (5) and (6) of reference I, the z,, here are not numerical 
values but functions of time while the subscript 7 refers not to the 
time of measurement but to the particular isotope employed. 

If the system being investigated is conservative 


N 


where K, is the amount present of the /th species. 
It follows that 


Viye*- Ci hp eayyee: @n 4 
V= |: EA og ea 3) 


Tine eTown Ky @i41,y e+ * enn 


Equations (3), (14), and (16) can therefore be used to determine 
the transfer rates per unit volume a, (0), the volumes V,(¢), and the 
transfer rates PAC) for an N compartment conservative system. 
Continuous access to all compartments is required throughout the 
course of the experiment in order to adequately determine z,,(¢). If 
the system under investigation is periodic or can be repeatedly 
prepared and a single isotope may be so introduced as to result in 
N independent initial conditions, solutions for a,, (¢) etc. can be 
obtained using only one isotope. In this event 2,,(¢) would repre- 

6As pointed out by Sheppard and Householder (Joc. cit.) only N=1 
radio-isotopes need be employed if a physiological substance is being 
studied. 
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sent the concentration of label in the ith compartment during the 
ith study. 

In certain special cases access to all compartments is not re- 
quired for some analytical purposes, and fewer isotopes may be 
needed. For a catenary system 


Git Fii41 Cini t Ceti + %,2-1%1-1,17 0, b= 1,2,3. (17) 
corresponding to equation 17 of reference I while the zth peripheral 
compartment of a mammillary system gives rise to the simple 


equation 
bi, + @;,(@,, — 24.) = 95 l=1 (18) 


In both cases the indicated a, ;(2) can be found if the time de- 
pendence of the appropriate concentrations is determinable. 


Il. Diffusion Problems. 


Since the concentration of labeled material is assumed uniform 
throughout a given compartment, it is almost a contradiction a 
priori to consider diffusion problems in terms of a finite number of 
compartments. In fact, in a rather definitive approach to the anal- 
ysis of the passage of a diffusible material through the walls of 
the vascular system, Sangren and Sheppard (1953) and Landahl 
(1954a) find little need to make use of compartment theory. None- 
theless, since one may conceive of the mammalian organism as 
consisting of a vascular compartment plus extravascular compart- 
ments and since this conception has validity for certain sub- 
stances such as inulin and iodinated albumin etc., it is only nat- 
ural to attempt to extend the area of application to materials which 
diffuse more readily through the walls of the vascular system and 
so gain some additional insight into the meaning of actual experi- 
mental results observed following intravenous injection of such 
diffusible substances. 

More precisely it would seem of interest to investigate the rela- 
tion between extravascular intracompartment concentration gradi- 
ents, intercompartment transfer rates, localization rates, and the 
observed experimental time dependence of the average concentra- 
tion in the one compartment (vascular compartment) in which mix- 
ing is usually rapid and complete. 

Consider the model illustrated in Figure 1 and assume that steady- 
state conditions prevail (i.€., 745, Tyo) Te1> 7.5 are constant). 
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If an isotope is introduced at point P, in compartment 1, the 
amount of material localized in 1. and 2. and exchanged pewveci 
1. and 2. will depend upon the concentrations along interfaces and 
not directly upon the average concentrations 21,'®% in 1. and 2. 


| toi B a 


oe 


fei lee 


fie 


. Se 


FIGURE 1. 


respectively.” Since, however, the label is present in tracer 
quantities only, the average concentrations z, must still obey 
linear equations of the form below. 


E,+ ) | 0,2, = 0. ear eae (19) 


pa 


The quantities a,, defined by equation (19) will not be constants 
but rather functions of the time following injection, the injection 
sites, and of course depend indirectly upon the diffusion and mix- 
ing characteristics within each compartment. Although very diffi- 
cult to determine, the a,, will generally be well-behaved functions 
of the time after injection 7, which approach the steady-state values 
az; The presence of intracompartment concentration gradients 
may therefore be thought to change a steady-state N compartment 
system into a kind of very complicated N compartment non-steady- 
state system. 
Examples of possible a,, are indicated below. 


7Determinations of the average concentration would of course be quite 
difficult in any compartment for which internal mixing is slow. Experi- 
mental values of the concentration will often be restricted to the vascular 


compartment. 
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FIGURE 2. 


From the form of the a;,; . it suggests that they can be very closely 
represented by: 


Gi, = 5; + » Aime *™, (20) 


m= 4 
where A,, are positive real constants. 


co 


+) Ame m7G,+ DAs e7 7B, = 0; (21) 
m =0 s=0 

+ D>) Agee E, + 3 A,o,@ 7 Z=0, (22) 
k=0 


where A io = 4:5) and X,, 
m=s=k=v. 

Expressed as second-order differential equations the coupled 
first-order system becomes: 


, A,, Ay, A, are not necessarily equal for 


Y, Fue eB, +) bge TE, +) ce M7 B, = 0; (28) 
a=0 B=o y=0 


Ye ase" hae F, + » tye AGTE fo i+ Sage e7Z, =0. (24) 


E€=Q 
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Equations (23) and (24) can be solved by series approximation 
to arbitrary accuracy for 7 sufficiently large, i.e., so that the series 
on the right side of equation (26) below converges. The method as 
well as the nature of the solutions will be illustrated by consider- 
ing a special case of equation (23): : 


Let x, + 1e4+ be] = + 1eyde 7] 2, =, 0, (25) 

and assume cece > A,,e7 aT , (26) 
w=0 

then 


Ss A? Ae hoT + e+ 5 As hala" 


@=0 
w 


a es ue] b Age" Ae C6) 


where 0 <A, cA5 Shy se SACS 


Letting 7 become large, it follows that 


[AZA, - aA, + cAgle *7 =0, (28) 
oe 
whence Ay = aoa tin ws (29) 


Since Z(7) is bounded and probably approaches equilibrium values 
, it seems reasonable that each of the roots of equation 
(29) corresponding to the two independent solutions of equation 
(25) are real and > 0.° Let 


Aw= K coe to Ls (30) 


2S 4% -——> 00 


It follows then that a recursion relation in the A, of equation 
(26) may be obtained by successively setting the coefficients of 
each exponential term in equation (27) equal to 0: 


_ (rons = DAww1 (31) 


2 
Ney GA +e 


A 


@ 


8In a detailed investigation of the relation between the nature of the 
roots and the transfer and localization parameters, this conclusion is 
verified for most steady-state systems (Schoenfeld, 1956). 
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By proper selection of a,b,c and d the numerator of (31) can be 
caused to vanish for w > some integer and the resulting complete 
solution of equation (25) may take the form of the sum of two in- 
finite series of exponential terms or the sum of one finite series 
of exponentials and one infinite series. 

In any event the recursion relation (31) and equation (26) show 
that the expansion of z, in terms of decaying exponentials (though 
rapidly convergent in this case for 7 > 9) can give rise to an infinite 
number of terms in steadystate systems with a finite number of 
physical or chemical compartments whenever intracompartment con- 
centration gradients exist. 

From the above it is clear that: 

(a) Tracer experiments with diffusible substances in steady-state 
systems will usually if not always give rise to experimental curves 
which can be broken down into or approximated by separate ex- 
ponential terms. 

(b) These separate experimental terms need not reflect the action 
of different physical processes or compartments but may only be 
the result of diffusion effects within the compartments. 

(c) As might have been expected, the effects of diffusion will be 
separable from those of the transfer and localization rates only if 
the time dependent exponential terms in a@,, decay more rapidly 


asl Sa : : 
than the e *°” terms where the A, are necessarily only functions 


of the equilibrium values a;; Otherwise non-negligible Pi 


—2 ‘ E = 
e 2” terms, etc., would be interspersed with e Xo” terms for all 


values of 7 and at present the theory of tracer analysis is not 
equipped to handle such a situation. 


III, Discussion. 


From the foregoing treatment it is apparent that tracer analysis 
of the kinetics of multicompartment biological systems must be 
approached with care. Since many different systems can give rise 
to much the same experimental results (especially when there is 
limited access to some of the compartments), the kinetic models 
obtained by interpreting tracer concentration curves should be 
investigated critically to see whether the assumptions concerning 
steady or non-steady-state conditions, intracompartment concentra- 
tions, and experimental accuracy are realistic. In Table I below, 
a classification of different types of commonly occurring systems 
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is introduced and appropriate analytical methods of study are 
indicated. 


TABLE I 

—————— 

Biological 
Classification Properties Occurrence Analytical Methods 
Steady-State ri¥i constant? Processes carried Multi-label analysis, 
Class 1 out under complete Multi-time analysis, 

steady-state con- For a general 
ditions. method of analysis 


of partially complete 
data, see Berman 
and Schoenfeld 


(1956). 
Steady-State rj; constant, Growth, pathologic Multi-label analysis 
Class 2°° V; variable conditions, adminis- (quasi-steady state 
tration of labeled treatment), Method 
non-physiological of Section I. 


substances which 
are localized. 


Non-steady r;j and/or Growth, pathological Method of Section I. 
state V; functions of conditions, response 
Class 3 time. to treatment, time 


delay systems. 


Pseudo non- raj and/or ef- Whenever compart- Diffusion equation 
steady state fective V; de- ment idealization analysis. 
Class 4 pendent upon breaks down and 
time after in- substantial intra- 
jection of compartment concen- 
label and or tration gradients are 
site of present. 
administration. 


It should be pointed out that the methods listed in Table I serve 
principally to evaluate parameters rather than to determine why 
they have the values they do. In order to examine this question, 
it will probably be necessary to study stimulus-response effects 
on non-tracer levels. As indicated by Danziger and Elmergreen 
(1954,1956), and others, the numerical complexity of such an in- 
vestigation suggests that modern digital and analog computers may 
play an important part in future developments. 


°See footnotes on page 88 of reference I. 
Steady-state Class 2 systems are so described because the metabolism 


of the labeled material remains unaltered even though the compartinent 
’ volumes are changing. 
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BOOK REVIEW 


Medical Research: A Midcentury Survey. Volume I: American Medical 
Research. Volume II: Unsolved Clinical Problems. Report of the 
American Foundation. 1955. 740 pp. Boston and Toronto: Little, 
Brown and Company, 


This two volume survey appears under the anonymous authorship of the 
American Foundation. Except for the foreword signed by Curtis Bok, 
outlining the general aims of the Foundation, and the introductory section 
on ‘‘Why and how this study was made’’ signed by Esther Everett Lape, 
who, together with a small staff, must evidently be credited with the 
compilation of the materials, the specific authorship of any of the essays 
comprising the work remains unknown. This is the more remarkable 
because, in spite of the great bulk of material on the level of factual 
reports comprising about three-quarters of Volume I, the work seems to 
have an author endowed with a keen understanding of the most fundamen- 
tal problems confronting the scientific community: the relation between 
pure and applied research; the problems arising from the conflicting needs 
of specialization and integration; the spectrum of organization of scien- 
tific effort, ranging from anarchy to bureaucracy; the search for talent; 
the tasks of communication; the relation between the financing of re- 
search and the control of its direction; and finally putting into proper 
perspective the great unsolved problems pertinent to medical research 
with which we are now confronted. 

Not only is the statement of all these problems lucid and forceful, 
but there is also a strong impression that the composite ‘‘author’’ comes 
to grips with these problems in a direct, positive way. In certain matters 
a definite stand is taken; in others, the pros and cons of a controversial 
issue are stated with admirable objectivity. The frequent citations serve 
to identify certain points of view as associated with certain personalities 
in American science. Nevertheless, the impression of a composite ef- 
fort is preserved throughout, leading to a hopeful conclusion that ‘‘Medi- 
cal Research’? in the United States has a mind, a voice, and an aware- 
ness of its role. 

Underlying the basic view of medical research as presented in the 
Survey is a clearly stated thesis, namely, that it is not sufficient to 
finance medical service to insure its availability, for clearly no amount 
of money will buy more medical service than can be got on the basis of 
available knowledge. The state of medical knowledge, in turn, depends 
on the quality of medical research, and the latter is limited by the state 
of the disciplines which feed it. These ‘‘disciplines’’ Comprise practi- 
cally all of science. First and foremost among the sciences relevant to 
medicine is certainly biology. But biological science, it is wisely 
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pointed out, is defined not so much by the departmental listings in uni- 
versity catalogues as by the kind of knowledge which is pertinent to the 
understanding of the life processes. This knowledge is derived certainly 
from chemistry and physics as much as from biology proper as it is con- 
ventionally understood. It is also derived from understanding phenomena 
usually assigned to social science, and also from mathematics. The 
importance of mathematics in biology, it is noted, is not limited to its 
role as the language in’ which practically all of physics and much of 
chemistry is written. Mathematics functions in biology, as in all science, 
as a focal point of theoretical organization. That this role of mathema- 
tics is clearly recognized in the survey is quite evident. There is a dis- 


novel approach to cell physiology through mathematical biophysics; of 
mathematical theories of the central nervous system. 

It is indeed gratifying to find not only a clear realization of the broad 
scientific basis which underlies medical knowledge but also the central 
position accorded to this realization. Thus a vantage point is obtained 
from which a survey of medical research and associated problems can be 
made in an extraordinarily meaningful way. What would otherwise be a 
bewildering compilation of reports, indices, and bibliographies becomes 
organized around well-discernible trends and central problems. 

The trends in medical research seem to be in the direction of integra- 
tion of concerted efforts by teams of workers from different traditionally 
academic ‘‘disciplines’’ toward an attack on what are felt to be central 
problems. An examination of these reveals a basic shift in the philos- 
ophy of medicine away from the ‘‘single cause’’ idea of disease which 
had permeated both the pre-scientific period and a considerable portion of 
the scientific era (the search for specific micro-organisms, later the 
search for specific diet deficiencies as the sole determinants of diseases). 
The abandonment of simplistic explanations in favor of an orientation to- 
ward multiple causes (presence of micro-organisms and immunity charac- 
teristics; hereditary and environmental factors) has not merely compli- 
cated the picture but has suggested a new focus of interest—the meta- 
bolic process, the derangement of which is the basis of disease. Thus a 
new unified view of pathology arises and necessitates a shift in contem- 
porary medical research from an exclusive preoccupation with the ab- 
normal to a systematic investigation of the normal functioning of the life 
process with which biology as a whole is concerned. 

To be sure, the sphere of biology must still be divided into areas of 
special interest, but these areas can no longer be determined by the 
traditional departmental divisions, such as botany, zoology, anatomy, 
and physiology. The re-grouping suggested by Weiss, for example, for 
administrative use within the National Research Council lists six major 
areas of biology: molecular, cellular, genetic, developmental, regulatory, 
and environmental. It is clear that this classification is a result of a 
recognition that the levels at which events occur and the fundamental 
processes of the living process ought to be the organizing principles of 
research. 

In line with these convictions, namely, that contributions to medical 
research ought to be determined by more profound considerations than the 
preservation of traditional departmental lines, the survey proceeds with 
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terse but content-packed accounts of the contributions to medical re- 
search from various non-biological sciences, mathematics, chemistry, 
and physics (including the development of atomic energy). It may be of 
special interest to the readers of this journal that the role of the various 
developments in bio-mathematics (mathematical biophysics, mathematical 
genetics, cybernetics) is amply recognized, 

Following an enlightening discussion of the relation of basic to clini- 
cal research, the survey raises the complex questions pertaining to the 
psychology and sociology of scientific research in general. It examines 
the current motivations for a research career, raises the problem of 
selecting personnel, and airs the great controversy of planned vs. in- 
dividually inspired research—the problem of ‘freedom vs. organization”? 
and of individual vs. project research. Both sides of this controversy are 
presented with admirable clarity and objectivity and the best is made of 
each case in the sense that really the strongest and the sincerest argu- 
ments are put forward. 

The important distinction between the ‘‘front line troops,’ i.e., re- 
search concerned with pushing the frontier of the unknown outward, and 
the “‘mopping up’’ and service echelons, i.e., research concerned with 
more routine investigations and applications, is brought out in support 
of the sober idea (only suggested) that different degrees of organization 
may be appropriate to each type of activity. 

There follows the sociology of research, the role of the foundations, of 
the national voluntary health organizations, of the pharmaceutical labora- 
tories, of the independently endowed research institutes, of the uni- 
versities and medical schools. The sociology of research is, of course, 
intimately linked with its economics, as reflected in the financing of re- 
search by philanthropic support, by industry, particularly the pharma- 
ceutical industry, and, of course, by the government whose role in the 
situation over just a decade or so has become of decisive importance. 

Finally the general social influences in relation to biological and 
medical research are discussed in relation to the status of the sciences 
and of medical practice and in relation to medical education. 

There follows the predominantly factual portion of the survey, the 
listing of the research agencies with compact resumes of their contribu- 
tions together with an account of the several financing arrangements. 
The recent research history and the present situation of some dozen 
major universities is so described. In addition, the work of the founda- 
tions, of non-teaching hospitals, of industry, and of government agencies 
is similarly treated. 

Aside from the first two chapters of Volume I, which raise crucial 
questions in the area of philosophy of science and sociology of knowl- 
edge, the scientifically oriented reader will be most attracted by the 
contents of Volume II. devoted to the discussion of the great unsolved 
problems of medicine. Following an introductory chapter on the ‘‘current 
metabolic concepts orienting research in biology and medicine,’’ the nine 
remaining chapters of Volume II are devoted to the following problems: 
cancer, infertility, arteriosclerosis, hypertension, the rheumatic syn- 
dromes, tuberculosis, viruses, alcoholism, and schizophrenia. Limitations 
of space preclude anything more than a listing of sub-topics of one or two 
of these areas which may give the reader an idea of the level of discus- 
sion. Cancer is discussed in the light of carcinogens, their action on 
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cells and tissues, their molecular structure, their role in mutation, hered- 
itary and somatic; in the light of hormones including the present status of 
sex hormone therapy, the role of the pituitary and adrenal hormones, and 
the problems related to the mechanism of action of hormones in metab- 
olism; in the light of enzymes and nucleoproteins, including the deter- 
mination of individual enzymes in research and diagnosis, etc; in the 
light..of the virus theory of cancer and its ramifications. 

The discussion of schizophrenia involves hereditary and age factors and 
concepts derived from emphasizing psychological, physiological, and social 
factors. The account of recent emphasis on physiological research con- 
tains an especially fascinating discussion of trends toward establishing 
psychological correlates of neurophysiological function and neuro- 
anatomical structure and arationale of attempts to apply neurophysiological 
principles to brain surgery. Biochemical psychological and cultural 
trends are discussed with equal thoroughness, 

This enumeration is, of course, no more than a paraphrase of a small 
part of the table of contents. This is to show that the compilation of an 
exhaustive table of contents of medical research would in itself be a 
notable achievement. The midcentury survey of medical research con- 
tributed by the American Foundation has exceeded this achievement 
many-fold by ingeniously combining in a single work a wealth of purely 
factual information on the status of medical research in the United States 
with a skillful and appealing formulation of an enlightened point of view. 
The book should be in possession of any who are interested or involved 
in medical research. 


ANATOL RAPOPORT 
Mental Health Research Institute 
University of Michigan 


THE SCIENTISTS’ COMMITTEE ON SECURITY 


The Scientists’ Committee on Security, Inc., an independent 
volunteer group formally organized on January 1, 1956, is now 
planning to expand its activities in the area of science and se- 
curity. The Committee seeks to fulfill several functions: 


1) To act as a clearinghouse for information and responsible sci- 
entific opinion on matters of information security and personnel 
security. The Committee attempts to answer enquiries from in- 
dividual scientists and others. Upon request, it will offer informal 
suggestions to help ensure that individuals with clearance diffi- 
culties receive the full protection of existing regulations, but it 
does not judge the merits of individual cases. 

2) To undertake, from time to time, special studies of particular, 
important issues of science and security and to prepare and make 
public appropriate reports. 

3) To cooperate with and to assist Government agencies in estab- 
lishing and maintaining realistic security programs which safe- 
guard both the long-range security of the United States and the 
traditional rights of its citizens. 

4) To stimulate constructive thought on questions of information 
security and personnel security. 

5) To foster a better popular understanding of the true relationship 
of science and security. 


Present members of the Scientists’ Committee on Security are: 
Ernest C. Pollard, Chairman; Ralph S. Brown, Jr., Vice-Chairman; 
John B. Phelps, Secretary; Hans G. Graetzer, Treasurer; Earle C. 
Fowler, Robert L. Gluckstern, Samuel A. Goudsmit, William A. 
Higinbotham, Franklin Hutchinson, G. Evelyn Hutchinson, William 
J. Knox, Henry L. Kraybill, Lester Lewis, Theodore K. Osgood, 
George F. Pieper, Waldo Rall, Julian M. Sturtevant, Hugh C. Wolfe. 

For its expenses, the Committee presently has a very small 
budget maintained entirely by voluntary contributions, mostly from 
working scientists. SCS members have written several articles in 
the area of science and security. Members have conferred informally 
with Government officials and have presented formal testimony 
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upon request. The Committee has submitted documentary material, 
with examples and illustrative personnel security cases, to Govern- 
ment groups on other occasions. It is expected that all these ac- 
tivities will be expanded. 

In recent months the Committee has learned of many instances 
in which clearance problems have affected the ability of scientists 
and engineers to find employment, commensurate with their ability 
and experience, outside the classified technical areas where clear- 
ance is understandably required. This situation—clearly contrary 
to the intent of security policies—is manifested in a variety of 
individual cases and is believed by the Committee to be serious. 
The Committee is now preparing a report on this special problem 
and will be grateful for any relevant information or examples. The 
Committee is also assisting the Commission on Government Se- 
curity and has an immediate need for detailed information on actual 
security cases involving scientists or engineers, for use in formu- 
lating and supporting recommendations for changes in personnel 
security regulations. All correspondence on individual cases is 
strictly confidential. The Committee earnestly solicits the com- 
ments and suggestions of scientists and others, and will appreciate 
views and information on security matters. Communications should 
be sent to the Scientists’ Committee on Security, Inc., 2153 Yale 
Station, New Haven, Connecticut. 
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Notice of Change of Title and Coverage 


INTERNATIONAL ABSTRACTS OF 
BIOLOGICAL SCIENCES 


Honorary Editor: C.C.N. Vass, M.Sc., Ph.D., M.B., Ch.B. 


The Council of Biological and Medical Abstracts Ltd. has decided to 
change the title of their journal British Abstracts of Medical Sciences to 


International Abstracts of Biological Sciences 


This decision results from two major developments. The first is an ex- 
tension of coverage to a wider range of biological research subjects pub- 
lished in journals throughout the world, and the second is that the Inter- 
national Abstracts of Biological Sciences, with the co-operation of the 
specialist editors selected by the Institute of Scientific Information of 
the Academy of Sciences of the U.S.S.R. will include, as from the January 
1957 issue, translations of the important Russian papers abstracted in 
the Referativny Zhurnal Biologii (Soviet Biological Abstracts) and 
Referativny Zhurnal Biologischekot Khimiit (Soviet Abstracts of Biological 
Chemistry), to be published simultaneously with their appearance in 
Russian. 

International Abstracts of Biological Sciences abstracts the world 
literature in the following fields: Anatomy, Animal Behaviour, Bio- 
chemistry, Biophysics, Cytology, Embryology, Endocrinology, Epidemiol- 
ogy, Experimental Biology, Genetics, Haematology, Histochemistry, 
Histology, Immunology, Microbiology, Nutrition, Odontology, Parasitology, 
Pathology, Pharmacology, Physical Anthropology, Physiology, Radiation 
Effects, Toxicology and Viruses. 

International Abstracts of Biological Sciences is published for the 
Council of Biological and Medical Abstracts Ltd. by Pergamon Press, 
London and New York, on a non profit-making basis. On the Council are 


represented: 


The Association for the Study of Animal Behaviour 
J. D. CARTHY, M.A., PH.D. 
The Biochemical Society 
R. K. CALLOW, M.A., D.PHIL., B.SC. 
R. A. ROBINSON, M.SC., LL.B., F.R.1.C. (Hon. Sec.). 
The Association of Clinical Pathologists 
A. G. SIGNY, M.B., B.S., M.R.C.S., L.R.C.P. 
The International Association for Dental Research 
A. E. W. MILES, F.D.S., L.R.C.P., M.R.C.S. 
The Society for Endocrinology 
I. W. ROWLANDS, M.SC., PH.D. 
The Society for Experimental Biology 
To be appointed. 
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The Nutrition Society 
G. H. BOURNE, D.PHIL.,.D:sG: 

The Pharmacological Society 
C.-A. KEELE, M.D., F.R.C.P. 

The Physiological Society 
H. P. GILDING, M.A., M.D. (Hon. Treasurer) 
W. D. M. PATON, M.A., D.M. 
C. C. N. VASS, M.SC., PH.D., M.B., CH.B. 
E. B. VERNEY, F.R.S., F.RSC.P; 


Since the Russian abstracts will have to be translated, the Council re- 
quires translators with knowledge of Anatomy, Animal Behaviour, Bio- 
chemistry, Endocrinology, Experimental Biology, Experimental Medicine, 
Microbiology, Odontology, Pathology, Pharmacology and Physiology, and 
would welcome offers of help. Terms of payment and other details may 
be obtained from the Honorary Editor: Dr. C. C. N. VASS, c/o Pergamon 
Press Ltd., 4 & 5, Fitzroy Square, London, W.1, or 122 East 55th Street, 


New York, 22, N.Y. 


